Gravitational Waves from Quasicircular Extreme Mass-Ratio Inspirals 
as Probes of Scalar- Tensor Theories 
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A stellar-mass compact object spiraling into a supermassive black hole, an extreme-mass-ratio 
inspiral, is one of the targets for future space-based gravitational-wave detectors. Such inspirals 
offer a unique opportunity to learn about astrophysics and test General Relativity in the strong- 
field. We here study whether scalar-tensor theories in asymptotically flat spacetimes can be further 
constrained with these inspirals. In the extreme-mass ratio limit, and assuming analyticity of the 
coupling functions entering the action, we show that all scalar-tensor theories universally reduce to 
massive or massless Brans-Dicke theory. We also show that in this limit, black holes do not emit 
dipolar radiation to all orders in post-Newtonian theory. For massless theories and quasi-circular 
orbits, we calculate the scalar energy flux in the test-particle, Teukolsky formalism to all orders in 
post-Newtonian theory and fit it to a high-order post-Newtonian expansion. We then derive the 
post-Newtonian corrections to the scalar-tensor modified, Fourier transform of the gravitational wave 
response function and map it to the parameterized post-Einsteinian framework. With the Teukolsky 
flux at hand, we use the effective-one-body framework adapted to extreme mass-ratio inspirals to 
calculate the scalar-tensor modifications to the gravitational waveform. We find that such corrections 
are smaller than those induced in the early inspiral of comparable-mass binaries, leading to projected 
bounds on the Brans-Dicke coupling parameter that are worse than current Solar System ones. This 
is because Brans-Dicke theory modifies the weak-field, leading to deviations in the energy fiux that 
are largest at small velocities. For massive theories, superradiance can produce resonances in the 
scalar energy fiux that can lead to quasi-circular floating orbits outside the innermost stable circular 
orbit and that last until the supermassive black hole loses enough mass and spin-angular momentum. 
If such floating orbits occur in the frequency band of a LISA-like mission, they would lead to a large 
dephasing (typically ~ 10^ rads) that would prevent detection of such modified inspirals using 
General Relativity templates. A detection that is consistent with General Relativity would then 
rule out the presence of floating resonances at frequencies lower than the lowest inspiral frequency 
observed, allowing for the strongest constraints yet on massive scalar tensor theories. 

PACS numbers: 04.80.Nn,04.30.-w,04.50.Kd 
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I. INTRODUCTION 

Gravitational waves hold the promise to allow for new 
and stringent tests of General Relativity (GR) in a pre- 
viously obscure regime: the strong-field region. This is 
the regime where gravitational fields are strong and the 
characteristic velocity of inspiraling bodies is a signifi- 
cant fraction of the speed of light. The inspiral, merger 
and ringdown of compact objects is the prototypical ex- 
ample of a gravitational wave source that emits in this 
region. We currently lack experimental verification of GR 
in the strong-field, as the most stringent tests we possess, 
thanks to the discovery of the double binary pulsar [l|, , 
only probe the weak field. For example, the ratio of the 
total mass to the orbital separation (the gravitational 
compactness), which controls the magnitude of the or- 
bital velocity, is of O(10~'^) <C 1 for the double binary 
pulsar while it is of 0(1) during binary coalescence. 

Much work has been devoted to studying whether grav- 
itational waves from ground-based and space-based de- 
tectors could be used to test GR. In most studies, one 
chooses a particular theory, or effect that one wishes 



to constrain, and then derives how gravitational waves 
would be deformed. Then, if a gravitational wave is de- 
tected that is consistent with GR, a constraint on gravi- 
tational wave deformations can be placed, which in turn 
are controlled by the magnitude of the coupling con- 
stants of the theory. Such a p rocedure has been applied 
to Brans-Dicke theory j3l-[l(j|. dynamical Ghern- Simons 
gravity [TT| - [l^ . modified quadratic gravity [l4|, phe - 
nomenological massive graviton propa gati on 10l.ll5l- 
[TtI . gravitational Lorentz-violation |18|. gravitational 
parity- violation [13, [l^ , violations of Local Position In 
variance 12011 and theories with extra-dimensions 21 



Most of these studies considered the early inspiral of 
compact objects within the post-Newtonian (PN) ap- 
proximation, where fields are expanded in small veloci- 
ties and weak-fields to high order in perturbation theory. 
Moreover, they typically considered only the leading- 
order (Newtonian) correction introduced by the given 
alternative theory and neglected all higher-order correc- 
tions. Such an approximation becomes particularly ill- 
suited in the late phase of inspiral, when the binary is 
close to the inner- most stable circular orbit (ISCO). This 
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then suggests that highly relativistic systems might allow 
us to place different, perhaps more stringent constraints 
on modified gravity theories. 

An example of such a highly relativistic system is an 
extreme mass-ratio inspiral (EMRI) , where a small com- 
pact object with mass mgco S (Ij 10)Mq spirals into a su- 
permassive black hole with mass Mmbh € (10^, 10^)AfQ. 
The EMRI phase of interest to future space-based grav- 
itational wave detection is the late inspiral, when the 
orbital separation is < 20Mmbh, the orbital period is 
(10^, 10'*) seconds and low-frequency gravitational waves 
are emitted in the (10~^, 10~^) Hz frequency band. The 
merger and ringdown EMRI phases are completely negli- 
gible, as they contribute insignificant amounts of signal- 
to- noise ratio to the signal [22|, . The EMRI inspiral 
alone can last tens to hundreds of years in the detector 
band, completely outlasting any realistic detection. In a 
single year of inspiral observation, millions of radians arc 
contained in the signal, which encodes rich information 
about the spacctime dynamics [24|. EMRIs arc there- 
fore exceptional probes of strong-field GR and one might 
expect them to be the best probes of strong-field GR 
modifications. 

In this paper, we study whether the detection of equa- 
torial and quasi-circular EMRIs allows for stronger con- 
straints on a large class of scalar-tensor theories than 
those currently in place due to Solar System observations 
or those projected with future gravitational wave obser- 
vations of comparable-mass binary inspirals. We choose 
to work with equatorial and quasi-circular EMRIs for two 
reasons. First, EMRI formation scenarios exist where 
the small compact object is in an equatorial and quasi- 
circular orbit by the time it enters the detector sensitivity 
band. One such model is that of [1^ , where a stellar-mass 
compact object is either created in an accretion disk sur- 
rounding a supermassive BH or is captured by the disk. 
The accretion disk is expected to be in the spin equa- 
torial plane inside a few hundred gravitational radii of 
the supermassive BH (26j . Thus, in this scenario by the 
time the small compact object reaches sufficiently small 
separations to produce detectable gravitational waves, it 
would be in an equatorial and quasi-circular orbit. Sec- 
ond, although such EMRIs are less likely than inclined 
and eccentric ones (since other EMRI formation channels 
are more favourable to generic EMRIs), the tests we will 
discuss below only require a single detection. Nonethe- 
less, the study presented here should lay the foundations 
for extensions to more generic orbits if needed. 

Scalar-tensor theories are prototypical models for non- 
minimally coupled scalar fields, which arc common ingre- 
dients in low-energy effective actions of quantum gravity 
models (see e.g. [231). Furthermore, scalar-tensor theo- 
ries can be shown to include f{R) gravity as a particular 
case (see e.g. [H, [2^ for some reviews), thus being rele- 
vant as potential models for dark energy. Although here 
we focus on asymptotically flat geometries, our results 
should also apply to f{R) theories, usually developed in 
asymptotically de-Sitter backgrounds, under the reason- 



able assumption that the cosmological background does 
not affect the local physics (see Ref. [s^l for a discussion 
on this issue). 

The most general, stationary, axisymmetric and vac- 
uum solution to generic scalar-tensor theories is the Kerr 
metric coupled to a constant scalar field. In fact, Sotiriou 
and Faraoni [s^l recently proved, under quite generic con- 
ditions, that such Kerr black holes are the endpoint of 
gravitational collapse in any scalar-tensor theory. This 
proof, however, fails for a cosmological, non-stationary 
(linear- in-time) scalar field [3l] |. which might lead to vi- 
olations of the no-hair theorem applied to scalar-tensor 
theories 

The Kerr metric will be adopted as the background 
spacctime upon which a small compact object inspirals. 
Although we mainly work in the Einstein frame, phys- 
ical observables arc obtained by mapping results to the 
physical Jordan frame. We show that, for a large class 
of scalar-tensor theories in an asymptotically fiat, sta- 
tionary background and expanding the fields about per- 
turbations proportional to the mass ratio, the equations 
governing the emission of gravitational and scalar waves 
only depend on two parameters: a coupling constant a 
and the mass of the scalar field ^s, both of which are 
uniquely determined by the specific scalar-tensor model. 
Thus, one can think of scalar-tensor theories as a contin- 
uous (Lie) space (manifold) of theories, where (a, /is) are 
coordinates. 

Brans-Dicke theory is contained in the subset of mass- 
less scalar-tensor theories, i.e. the subset of scalar-tensor 
theories with /is = 0. This theory not only requires a 
massless scalar field, but also that a be related to Wqd 
by a specific relation (cf. Eq. ([^^ below). Brans-Dicke 
theory is then parameterized by a single quantity, cjbd, 
and can thus be thought of as a line in the (a, /is) space 
of all scalar-tensor theories. 



Executive Summary of Results 

We study EMRIs in generic scalar tensor theories. By 
assuming analyticity of the coupling functions entering 
the action, we show that any of such scalar-tensor theo- 
ries in the EMRI limit reduces identically and universally 
to massless or massive Brans-Dicke theory. To next-order 
in the mass ratio, this degeneracy in theory space breaks 
down. Therefore, any bound on Brans-Dicke theory de- 
rived from EMRI calculations automatically constrains 
the linear-in-mass-ratio expansion of all scalar-tensor the- 
ories. We also show that in the EMRI limit, black hole 
binaries in any scalar-tensor theory do not lead to dipo- 
lar radiation to all orders in FN theory, provided the 
background scalar field is constant at spatial infinity. 

We then study massless scalar-tensor theories. First, 
we numerically calculate the scalar energy fiux carried to 
spatial infinity in the point-particle approximation within 
the Tcukolsky formalism and assuming quasi-circular or- 
bits. This calculation of the fiux is fully relativistic, al- 
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beit valid only to leading order in the mass ratio, but 
accounting for all orders in PN theory. We fit this flux to 
a 3PN and a 3.5PN expansion and calculate the Fourier 
transform of the gravitational wave response function in 
the stationary-phase approximation beyond Newtonian 
order. We find that the PN corrections to the mass- 
less scalar-tensor modifications to the Fourier phase can 
be mapped to the recently proposed parameterized post- 
Einsteinian (ppE) framework [22l.[23| in a straightforward 
manner, which we explicitly provide. 

We then proceed to determine whether these PN cor- 
rections to the modified gravitational wave allow us 
to constrain the Brans-Dicke coupling parameter better 
than with comparable-mass binary inspirals. Recall that 
due to the equivalence discussed earlier, the bounds de- 
rived here also apply to all massless scalar-tensor theories 
to leading-order in the mass ratio. With the Teukolsky 
flux at hand, we employ the effective-one-body (BOB) 
framework [sj, HJ] , recently adapted to EMRIs [s^ [l^ , 
to evolve inspiral with and without the Brans-Dicke cor- 
rection. We then compare their associated gravitational 
wave phases, after minimizing the difference over an ar- 
bitrary time and phase offset. Contrary to our expecta- 
tion, we find that EMRIs will not be able to constrain 
Brans-Dicke theory beyond current Solar System con- 
straints. In fact, the projected bounds we obtain are 
worse than those found by considering gravitational- wave 
detection of comparable-mass binary inspirals with only 
the leading-order Brans-Dicke correction 0, 0, l6l-[Tol| . 

This seemingly surprising result is in fact perfectly 
reasonable, once one realizes that scalar-tensor theories 
are not necessarily strong-curvature modifications of GR. 
At the level of the action, the scalar-field modification 
acts on the Ricci scalar in the Jordan frame and it 
does not introduce higher curvature corrections. This 
is consistent with the fact that scalar-tensor theories in- 
clude f{R) gravity, which is usually studied as an in- 
frared, rather than an ultraviolet (or strong- ficld/strong- 
curvaturc) correction to GR. In terms of the energy fiux, 
the weak-field nature of massless scalar-tensor theories is 
evidenced by the scalar flux becoming less and less im- 
portant than the GR one as the small compact object ap- 
proaches the horizon. This is why the modified theory in- 
troduces a pre-Ncwtonian correction to the gravitational 
wave phase (a — IPN effect relative to GR), instead of a 
posi-Newtonian one. 

With the massless case under control, we then proceed 
to consider massive scalar-tensor theories, where now GR 
deformations depend both on the mass of the scalar field 
fis and the scalar-tensor coupling parameter a. Depend- 
ing on the values of {a, fig), prograde EMRI orbits can 
experience superradiance and produce resonances in the 
scalar fiux that can momentarily counteract the gravita- 
tional wave flux, leading to floating orbits [s^ , i.e. orbits 
where the inspiral greatly slows down at a given floating 
radius. Figure [1] shows the scparatrix (red dashed line) 
between floating and non-floating orbits for an EMRI 
with total mass M ^ lO^MfT, in massive Brans-Dicke 




FIG. 1: Parameter space of massive Brans-Dicke theory for 
a quasi-circular EMRI with typical neutron star sensitivity 
Ssco = 0.188, supermassive black hole spin aMBn/AfMBH = 0.9 
and mass Mmbh = 10"" Mq. The red, dashed curve is the 
boundary of the region in (lj-bu, /is) space that separates float- 
ing and non-floating resonances (cf. Eq. ((69])). The area below 
the black curve shows the region that is ruled out by Solar 
System experiments [s^, H^- The area between the green, 
dot-dashed curves shows the region where a resonance would 
occur inside the classic LISA sensitivity band [(10~*, 10~^) 
Hz]. Due to a large dephasing introduced by floating res- 
onances, any detection of an EMRI in a LISA-like mission 
would rule out the (open) region delimited by ABC. The dot- 
dashed blue line shows the values of (a;BD,/is) below which 
non-floating resonances lead to dephasings larger than 1 rad 
(see Sec. IVII"C)) . The ISCO frequency, fiisco, is shown by a 



vertical dotted line. 



The curves ^^^bd ^-nd ^^bd terminate at 
0.22/Mmbh. 



theory. As indicated in this figure, floating orbits can 
exist only in the region above and to the right of the 
red-dashed line. This curve (and all others) terminate 
when the small compact object approaches the ISCO, ap- 
proximately when /is ^ f^iscoi where r^isco is the ISCO 
frequency (vertical dotted line in Fig. [I}. 

We have chosen to show the region of allowed floating 
orbits in (ojbdi/^s) coordinates instead of (a,fis), thus 
restricting attention to massive Brans-Dicke theory via 
Eq. (p9)) . so as to be able to compare directly with So- 
lar System experiments [13, The area under the 
solid black line in Fig. [T] denotes the region that is cur- 
rently ruled out by such experiments [33, . Observe 
that as /is — >■ 0, the Solar System constraint becomes 
independent of /is and approaches uj^d > 4 x 10"*, as 
expected. For completeness, wc also show the region 
(between green, dot-dashed vertical lines) for which the 
floating orbits would lead to GWs inside the classic-LISA 
sensitivity band [(10"'', 10~^) Hz], since the resonant fre- 



quency is linearly proportional to /is (see Ref. [41| for 
a recent review on eLISA-NGO). From this, one can 
easily observe that there is a large region in (wbdjMs) 
space that is not ruled out by Solar System experiments, 
yet where floating orbits should be observed with LISA 
like instruments, ie. the region between the greened dot- 
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dashed curves, above the black sohd hne and below the 
red dashed line. As we explain below, the modifications 
to an EMRI GW signal due to hitting a fioating orbit 
are so large that an observation of a signal consistent 
with GR would rule out the region below the red dashed 
line, thus improving Solar System constraints by several 
orders of magnitude. 

Floating orbits are not a new invention in the context 
of EMRIs and GR. Press and Teukolsky [H, [H], after 
a suggestion of Misner's, were the first to study float- 
ing orbits within GR. The main idea is that for certain 
orbital configurations, a small object orbiting a super- 
massive black hole can extract angular momentum from 
the latter. That is, the gravitational waves that impinge 
on the supcrmassive black hole horizon are not absorbed, 
but instead refiected with greater energy and angular mo- 
mentum that they possessed initially. If such extraction 
is sufficiently large to balance the energy lost by the sys- 
tem due to gravitational waves traveling out to spatial in- 
finity, the small object might momentarily experience no 
radiation-reaction force and stall in its orbit. Press and 
Teukolsky found that, although such fioating orbits are 
not impossible in GR, for them to exist outside the ISCO 
of a prograde orbit, one requires an essentially maximally 
spinning supcrmassive black hole, with spins exceeding 
the Thorne limit |^max| = 0.9987^^^^^ In mas- 

sive scalar-tensor theory, more energy can be extracted 
due to the activation of the scalar energy flux, thus allow- 
ing for floating orbits outside the ISCO for supcrmassive 
black holes with smaller spins. 

Floating orbits, however, are somewhat of a misnomer, 
as an EMRI that attempts to cross such a floating radius 
does not completely stall or stagnate as the name implies. 
Instead, as an EMRI attempts to cross a floating radius, 
its evolution is greatly slowed down, as it is now driven 
by the supcrmassive black hole's energy and spin angular 
momentum loss. In other words, the radiation-reaction 
force or the energy flux are weakened by an additional 
factor of mass-ratio v. Thus, while to leading (first) order 
the energy flux scales with i/^, during fioating it scales 
with z/^, i.e. it becomes a second-order in the mass ratio 
effect. 

The time it takes the EMRI to cross a fioating res- 
onance is determined by the supcrmassive black hole's 
shedding efficiency [of 0(10^) yrs for a typical EMRI]. 
This timescale can easily exceed a Hubble time for EM- 
RIs orbiting at separations greater than 326AfMBH- While 
the EMRI is slowed down, it will produce gravitational 
waves very different from those in GR, leading to a 
large dephasing, i.e. a phase difference relative to the 
GR phase, of roughly 10^ rads over a one year obser- 
vation for phenomenologically viable {a, fig) parameters. 
Therefore, if massive, scalar-tensor theories are the cor- 
rect representation of Nature, such a large dephasing 
would prevent detection of EMRIs in a GR template- 
based, matched-filtering gravitational wave search. 

On the other hand, the detection of an EMRI grav- 
itational wave that is consistent with GR can then be 



used to constrain massive scalar-tensor theories. This is 
because such a detection would imply that a resonance 
was not present in the range of frequencies sampled by 
the EMRI. If a resonance were present at a frequency 
lower than the lowest one detected, the EMRI would 
have hanged at the resonance frequency for a very long 
time, evolving very slowly, and never entering the detec- 
tor's sensitivity band. The requirement that the resonant 
frequency be lower than the lowest EMRI frequency de- 
tected would impose unprecedented constraints on mas- 
sive scalar tensor theory: the (open) region delimited by 
ABC in Fig. [1] would be ruled out. Observe that such 
a constraint is much stronger than current Solar-System 
ones (solid black line in Fig. [T]). Therefore, even though 
future space-based gravitational wave detectors will not 
be able to constrain massless scalar-tensor theories be- 
yond current Solar System levels [stI . l40j . they will be 
able to severely bound a large sector of massive scalar 
tensor theories that is today still phenomenologically vi- 
able. 

The remainder of this paper presents the details of 
this calculation. Section lU reviews the basics of scalar- 
tensor theories applicable to EMRIs. Section IIIII de- 
scribes the test-particle approximation and the Teukolsky 
framework used to model EMRIs in scalar-tensor theo- 
ries. Section HVl calculates the energy fiux in the Teukol- 
sky formalism for massless scalar fields, fits the massless 
fiux to a PN expansion and computes the Fourier trans- 
form of the response function in the stationary-phase ap- 
proximation. Section |V] repeats the analysis of Sec. IIVI 
but for massive scalar fields. Section I VII describes how 
to model EMRIs using the EOB-EMRI framework both 
in GR and in massless scalar-tensor theories, and then 
calculates the dephasing between a GR waveform and 
a massless scalar-tensor theory one, allowing us to de- 
rive a projected bound. Section IVIII derives analytical 
estimates on how well massive scalar-tensor theories will 
be constrained with future gravitational wave detectors. 
Section IVlIII concludes and points to future research. 

Throughout the rest of this paper we follow mainly the 
conventions of Misner, Thorne and Wheeler [4^. We use 
Greek letters to denote spacetime indices, Latin ones at 
the middle of the alphabet i, j, . . . stand for spatial indices 
only, and Latin ones at the beginning of the alphabet 
a, 6, . . . range over waveform parameters. We also use 
geometric units with G = c = H = 1, except in Section HIl 
where we restore G for clarity. Here G is the gravitational 
constant as measured by Cavendish-like experiments, i.e. 
as measured by a distant observer in the Jordan frame. 



II. SCALAR-TENSOR THEORIES FOR 
EXTREME MASS-RATIO INSPIRALS 

In this section, we review the basics of scalar-tensor 
theories as they concern this paper. For a detailed review 
refer to [l^, S] and references therein. Consider then 
the action of a generic scalar tensor theory in the Jordan 
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frame: 



5, 



(J) 



1 

' 16^ 
-2C7(^) 



d^x^[F{cf,)R - Z{4>)g^^''d^<t>d,(l) 

+ <5'raat(*m;5Mi^) , (1) 



where F{(f)), Z{(j)) and [/(</>) are functionals of the scalar 
field (p that modify the Einstein-Hilbert term, the mini- 
mal kinetic term and the scalar field potential. The quan- 
tity <S'mat(^'m; ff/^y) dcnotcs tlic action for other matter 
degrees of freedom. As usual, R is the Ricci scalar and g 
is the determinant of the physical metric . 

In the following, we employ a test-particle approxi- 
mation, i.e. we assume the symmetric mass ratio i> = 
Tosco AfMBH/(™sco + Mmbh)^ < 1, where msco is the 
physical mass of a small compact object, while Mmbh 
is the physical mass of the supermassive (background) 
black hole. Thus, we consider the matter action for a 
point-particle moving on a fixed background 



Srr 



dr m{(j}) 



(2) 



where r is the proper time along the trajectory. Scalar- 
tensor theories in the Jordan frame promote the mass of 
small object to a 0-depcndent quantity m = m{4>) (47l . 
l48j . This allows for possible effects on the structure of 
the c omp act object, and on its motion, due to the scalar 
field glliilll. 

Let us now transform this action to the Einstein frame 
through the following conformal transformation (FH 



1 



3i^'(0)2 , 1Z(0) 



4 i^(0)2 2F(0) 



1/2 



The action in Einstein frame then becomes 



(3) 

(4) 
(5) 
(6) 



'(_E) / V a 1 

- / dT<'^^ A{<i>)m{(l)) 



tensor associated with the matter action, and 



where we have used dr = A(^)dT^-^'> . 
modified field equations are [50[ 



□ ^ ^ , 

167r 9$ yil^ 



9$ 



The associated 



(8) 



where denotes the location of the test-particle, 
j^f^AE) ^ 2(-5(^))-i/2,5S',„at/^5^<f^ is the stress-energy 



1 



is the scalar field energy-momentum tensor. Notice that 
the field equations depend only on three generic func- 
tions, and to(0). 

Let us now consider scalar perturbations around a 
constant background scalar field, and around a 

metric which is solution of Einstein's equations in an 
asymptotically-flat spacetime. Hence, we expand about 
<j) _ $(0) = ^ 1 assuming^ a general analytical behav- 
ior of the potentials around $ ^ $(0); 



n=0 

^ J2 ^« ~ 

?i=0 

m(<i>) 



n=0 



(9) 
(10) 
(11) 



where (Vn, An, nin) are constant coefficients independent 
of $ and m = to[0($)] = to($). Expanding the modified 
field equations in this way, we find 



□ 



where 



1^ 



^ /If 

167r 



1 



mi 
Too 



[i/o + yi(^] (12) 

Ti'(^). (13) 



(14) 

is the trace of the stress-energy tensor for a test particle 
in the Einstein frame. We have here dropped terms of 
0{lp v) and the sub or superscripts or 1 refer to quanti- 
ties that are truncated at zero or first order, respectively. 
Asymptotically flat backgrounds require Vo = = Vi, or 
equivalently Uq = Q = Ui. 

Our flnal equations then read 



167r J 
(7) 



GUE) 



where we have defined 

- ^1 ^ 
a = — 4- 

^0 



mi 
Too 



2_ ^2 



(15) 
(16) 

(17) 



^ By assuming V('I>), j4(<E>) and m(<i>) to be analytic functions of 
<I> we exclude some viable f{R) cosmological models which are 
equivalent to scalar-tensor theories with a divergent potential, 
— >• <E>(0)) — >■ oo (see e.g. Ref. [1^). However, in the context 
of strong-field corrections to GR, the analyticity assumption is 
general enough to encompass a large (in fact, infinitely dimen- 
sional) class of viable theories in asymptotically flat spacetime. 
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which are both constant. Notice that this is nothing but 
the Einstein equations for a particle with renormalized 
mass Tosco = Aqitio. Thus, any measurement of the 
mass of the particle is in reality a measurement of the 
renormalized mass TOsco j instead of the bare mass toq in 
the Jordan frame. 

Given this, the gravitational constant as measured by 
Cavendish type experiments (i.e. by observers at infinity 
in the Jordan frame) reads [49| 



G 



1 (A'im'' 



whereas the "sensitivity" of the compact object is 



d log TO 

dlogG 

TTTOiG 

AiniQ 



G($) m'($) 
to($) G'($) 



{A2+27TAoy 



(19) 



In the next sections, following the convention of previ- 
ous works [53, [12 , we shall use units such that G = 1 
which, by Eq. (HH), implies A§ = 1 - A^/(47r). This then 
allows us to simplify the expressions for the sensitivities 
in Eq. ([T9| and the a parameter in Eq. ([TTl) . 

If the small compact object is a black hole, then its 

Using 



mass TObh in the Einstein frame is constant [5 
the mapping to the Jordan frame, TObh 



mil^/A{^), 



one then finds 



TObh (^0 - Aiip) 



TOl 

Too 



and the black hole sensitivity reads 



\{A^ 







A2AQ 
27r 



4, (20) 



(21) 



Although Sbh generically depends on the specific scalar- 
tensor theory and can be different from its value in Brans- 
Dicke theory [s^ (sbh = 1/2, as discussed below) if the 
small compact object is a black hole we get, by Eqs. pT)) 
and (PO)) . a = 0. Thus, for any (analytic) scalar-tensor 
theory, scalar perturbations decouple from the matter 
sector and the energy and angular momentum emission is 
simply governed by Eq. (|15p , which is the same as in GR. 
Notice that this result is valid to all orders in PN theory, 
but only to leading order in the mass ratio. However, 
this result has been recently confirmed for equal-mass 
binary black hole mergers with full numerical relativity 
simulations [53 |. 

After choosing G = 1, the gravitational and scalar 
wave emission in the test-particle approximation only 
depends on two parameters, a and ^s, which are fixed 
by the specific scalar-tensor theory under consideration. 
Physically, a is related to the coupling function between 
the scalar field and the gravity sector, whereas is the 



mass of the scalar field. The latter is related to the pa- 
rameter of the coupling functions in the Jordan frame 
via 



.42^'($(o)) 



2V27r 



(22) 



where TOg is defined by U{(j)) = to^((/)— (/)o)^/2. Note that, 
due to the non-minimal kinetic term Z(0), the physical 
mass of the scalar field is ^s, and not TOs (see e.g. Ref. [sl] 
for a related discussion in massive Brans-Dicke theory). 

The modified field equations of scalar-tensor theories 
in the Einstein frame ([T5|) - ([TE1) are equivalent to those 
arising in GR coupled to a free scalar field if one expands 
in scalar perturbations and in the mass ratio. Thus, the 
gravitational metric perturbation only has the usual two 
propagating degrees of freedom. This is in contrast to 
modified field equations in Jordan frame, where the met- 
ric perturbation has three degrees of freedom: the two 
usual Einstein ones plus a breathing mode [52 1. 

In the next sections, we will study EMRIs directly in 
the Einstein frame and, once a solution is found, we will 
map the observables back to the Jordan frame. In the 
latter, the gravitational field at infinity reads 

5,. = A(a>)V5 = A($)2(^(f) + /^(f)) 



Ali^jlf^^ + hif) + 2A,Ao^4^J 



9^1 



(23) 



where we have neglected terms of Oi^pv) and quantities 
without superscript are in the Jordan frame. When com- 
puting waveforms, we are interested in the transverse- 
traceless part of g^i, only, hence the only relevant term 
in Eq. 



"'fJ.U 



2, TT,(£) 



Aih 



The quantity Aq is a constant that modifies the am- 
plitude of the transverse-traceless, metric perturbation, 
and thus, the gravitational fluxes of energy and angular 
momentum and the gravitational wave phase evolution. 
As discussed below, in massless Brans-Dicke theory, this 
is very close to unity: Aq = [(2a;BD + 3)/(2wbd + 4)]^/^ ^ 
1 — (4a;BD)~"'^, where Wbd is the Brans-Dicke coupling pa- 
rameter and in the last line we expanded about Wbd = 00 
since Wbd > 4 x 10^ to agree with Solar System experi- 
ments. Even in massive Brans-Dicke theory, where Wbd 
can evade Solar system constraints for a sufficiently large 
mass ^5, Aq is still close to unity: if cjbd = 1, then 
Ao = v%/3 ~ 0.745. 

Regardless of its magnitude, the quantity ^0 has the 
effect of renormalizing the overall amplitude of the grav- 
itational metric perturbation. This occurs because in 
Eq. (|T8|) we have implicitly chosen G^^' = 1 in the Ein- 
stein frame. We could have, however, chosen G*-^' in 
such a way as to force Aq = 1 exactly. In any case, for 
an EMRI observation, such an overall {G^^^Aq)-^ factor 
cannot be separately measured; it is degenerate with, for 
example, the luminosity distance. For this reason, we 
will ignore this overall factor of Aq in remaining sections 
of this paper. 
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Massless Brans-Dicke Limit 

Brans-Dicke theory can be recovered from the field 
equations ((T5|) - P^ (i.e. from the field equations in the 
test-particle approximation) in the massless, /i^ = 0, case 
with the choices 



F(0) = 0, Z{4>) 



U{^)=0, (24) 



where Wbd is a constant. This then leads to the conformal 
transformation ISOi 



Am 



9tii^ , 



1 

rj.(E) 



■ ln( 



-T„, 



(25) 



where /3 = yJVo-Kj (2a;BD + 3). In this case, the perturba- 
tion quantities Ai read 



/3 



^2 



and Eq. ([T8|) reduces to 



G 



00 2a;BD + 3 



(26) 



(27) 



Note that, if G = 1 units are used, the equation above 
defines 0o and Aq in terms of : 



2Wbd 



2a;BD + 3 



= An 



(28) 



Finally, the field equations in Einstein frame for an 
asymptotically-flat spacetime are given by Eqs. (jisp 
and ([161) with 



167rG 



2wbd + 3 

where the sensitivity of Eq. ([TO]) reduces to 

_ 1 mi 

Ssco = -T. , 



Ms = . 



(29) 



(30) 



which agrees with the results of Rcfs. [46!,l50| . If the small 
compact object is a black hole we get, by Eq. (|2T|) . Sbh ~ 
1/2. Finally, notice that Eq. ^ reduces to Eq. (2.11) 
in Ref. 

Brans-Dicke theory has already been strongly con- 
strained by Solar System experiments. GR is recovered 
in the Wbd ^ oo limit and the tracking of the Cassini 
spacecraft has lead to the constraints ojbd > 4 x 10* [40| . 



Massive Brans-Dicke Limit 

Massive Brans-Dicke theory can be obtained with the 
same choices of coupling functions F{(f>) and Z((j)) as in 



Eq. ([24|) . which then leads to the same choice of a as in 
Eq. (P^ . Massive theories, however, require the addition 
of a non- vanishing potential 



C/(0) 



2 



and hence, by Eq. 



^2 = and 



fJ-s 



2tJBD + 3 



(31) 



(32) 



Therefore, massive Brans-Dicke theory is described by 
two-parameters only, (ojbdi/^s)- 

Before proceeding, notice that scalar-tensor theories 
comprise a class of models, of which Brans-Dicke theory 
is in principle a single member. However, to linear or- 
der in the mass ratio and assuming analyticity, we have 
shown that all scalar tensor theories without (with) a po- 
tential reduce to massless (massive) Brans-Dicke. There- 
fore, any bound on the Brans-Dicke coupling parame- 
ter Wbd (or the Brans-Dicke mass ^s) is automatically a 
constraint on the coupling parameters of massless scalar- 
tensor theories (or the scalar mass) to linear order in the 
mass ratio. To next order in the mass ratio, this equiva- 
lence breaks down. 



III. EMRIS IN THE TEUKOLSKY 
FRAMEWORK 

In this section, we describe the Teukolsky framework, 
which we use to study EMRIs in the test particle ap- 
proximation. Hereafter, we restore the usual G = 1 
units. Consider then quasi-circular inspirals in the ex- 
treme mass-ratio limit: a small compact object with 
mass msco and no spin orbiting around a supermassive 
black hole with mass Mmbh, spin angular momentum 
IiS^bhI = Qmbh-^mbh and dimensionless Kerr spin param- 
eter gMBH = aMBn/MMBH- Wc rccah that mgco = Aqtuq, 
where mo is the bare mass. The total mass of the sys- 
tem is then M = msco + .^^mbh, the reduced mass is 
Msco = Tosco -^mbh/.^ ^ 1 and the symmetric mass- 
ratio IB u = ^sco/M ^1/4. We consider quasi-circular, 
equatorial EMRIs for simplicity. The spin of the small 
compact object can be neglected, as its effect on the evo- 
lution will be second-order in the mass ratio v 1. 
Finally, in the extreme mass-ratio limit, possible tidal 
effects may also be neglected [ssf . 

With this at hand, let us consider the Teukolsky for- 
malism [56j . In the test-particle limit, we can approxi- 
mate the EMRI trajectories as geodesies of a test particle 
with mass is a background of mass M to leading order 
in the mass-ratio. We can then study the gravitational 
and scalar waves emitted and the energy-momentum car- 
ried away by solving the first-order perturbation of the 
field equations as a function of the given geodesic. No- 
tice that the scalar-tensor modified field equations 
(1161) are very similar to the Einstein equations with a 
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propagatingscalar field, a system originally studied by 
Dctweiler j57| . 

Any given geodesic is extremely sensitive to the back- 
ground upon which it evolves, where we here choose the 
latter to be the Kerr metric. Sotiriou and Faraoni [s^l 
recently showed that the most general, stationary, ax- 
isymmetric, vacuum spacetime that is also the endpoint 
of gravitational collapse in scalar-tensor theories is the 
Kerr metric. This does not apply to non-stationary 
back grou nds, with time-dependent hair of the form ip ~ 
Mc t |3l| . where fic is a constant. In that case, ip = t^c 
is assumed to be of cosmological origin, and thus much 
smaller than any frequency in the EMRI problem. 

Therefore, without loss of generality, we linearize 
Eq. (|15p . the equations for the gravitational perturba- 
tions, about the Kerr background [in Boyer-Lindquist 
coordinates (t, r, 9, (/))] to first order in the mass ratio. 
One can decouple these equations in favor of a differen- 
tial equation for the Newman- Penrose scalar ^4, which 
we can harmonically decompose via 



tp4 



1 



(r - iA/MBH^MBH cos6')4 

/OO ^ 
d^^i?,™„(r) Si^UO)e'^""^-^'^ , (33) 

-°° Lni 

where the radial function Rimuj rnust satisfy the radial 
equation 



while the angular function Simu must satisfy the angular 
equation 



1 d 

+ 4 — 4m cos 



(^sinO-^SirnS^ - (^a^,j3jjW^ sin^ 6* + 4aMBH'^ cos 6* 
2aMBH(^ni - 2 - XjSim^ = (35) 



sm 



where {l,m) are harmonic indices, A is the eigenvalue of 
Eq. dSSl) and 



A = - 2MMBHr + M^^^ql^^, 



dr V dr 



K 



(36) 



The source term Timu can be derived from the stress- 
energy tensor of test-particles 



The equation for the evolution of the scalar perturba- 
tion in Eq. can be solved by decomposing the latter 
via 



ip{t, r, 6*, 0) = ^ j duje 



i7nd)—iujt 



(38) 

where again {l,m) are harmonic indices, Si„i{9) are 
spheroidal harmonics and Ximuiir) is a radial function. 
With this parameterization, the scalar field evolution 
equation becomes 



d^ 

drl 



A 



('^2 I q2 -j3/: 

V ' ""MBH/ 



(39) 



where dr/dr^ = A/ (?'^ + 0^^^) and the effective potential 
reads (see e.g. Ref. (50| ) 



Vs 



+aMBH(''^ - 4A/MBHr + aLn)] 



A/1 



(40) 



and A = Ai„, + al.^^u'^ - al^^^nl - 



2aM-B¥imuj is found from the angular eigenfunction 



1 d ( .^^gdS_ 
sh[ed0 \ 'de 



+ Ai, 



5 = 0. 



The scalar source function Ti„iuj is different from that 
which sources gravitational perturbations, Timu, and it 
is given by 



T^hn (■^) ^('' ~ ^sco) msco<5(TOSlorb - i^) ■ 



(41) 

where the overhead star stands for complex conjugation. 
We recall that a is given by Eq. (|29p in the Brans-Dicke 
case. 

Let us now consider two independent solutions Xi^^ 
and X^^^ to the homogeneous version of Eq. ([M)) . satis- 
fying the following boundary conditions: 



as r — > 00, , 



where fc+ = w — mfl^, f2+ — gMBH/(2r-|_) and koo = 
^yu}'^ — fxl- The quantity r+ is the largest root of A = 
0, namely r± = Mmbh ± -Mmbh (l ^ 'Zmbh); which is the 
location of the event horizon. Given this, the fluxes of 
scalar energy at the horizon and at infinity are [s^ 



-gu' 



-S{r-rsco)S{0-TT/2)S{^~n^,J) , (37) 



where rsco is the location of the small compact object, 
u'^ is the small compact object's four- velocity and fl^rb is 
its orbital frequency. The source term and the details of 
how to solve Eas. in terms of Green's functions 

are given in j57| . so we omit them here. 



Ini 



(42) 



^ Note that A is different from the usual A defined, e.g. by De- 
tweiler [5Sj and Teukolsky [59j (in their notation, A = A -|- 
'^mbh/^s)- 
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where we have defined 



X 



Irnuj 



(^sco) 



(43) 



and where W is the Wronskian of the two linearly inde- 
pendent homogeneous solutions, W = Xi^^dX^^/di\ — 

Through this formalism, we can then obtain the grav- 
itational wave and scalar field energy flux to all orders 
in PN theory. This is because the Teukolsky approach 
employed here only linearizes the field equations in the 
mass-ratio, while retaining the full relativistic nature of 
the orbits. This approach, however, is only valid for 
closed geodesies, i.e. provided the implicit averaging car- 
ried out through the harmonic decomposition is valid. 
This is sometimes referred to as the adiabatic approxi- 
mation, where one assumes the worldline of the small 
compact object can be approximated as a sequence of 
osculating geodesies. Such an approximation, of course, 
breaks down when the small compact object enters the 
ISCO, which is where we stop all our evolutions. Fortu- 
nately, the missing plunge and ringdown phases are com- 
pletely negligible for data analysis purposes given EM- 
RIs, since these two phases contribute negligible to the 
total signal-to-noise ratio. 



IV. MASSLESS SCALAR ENERGY FLUX 

Let us first concentrate on massless scalar-tensor the- 
ories and compute the emitted scalar energy flux. Be- 
fore solving the equations of Sec. IIIII numerically, let us 
describe certain analytical solution in the PN approxi- 
mation, V <^ 1 where v is the binary's relative orbital 
velocity. In this limit, the scalar fiux carried out to spa- 
tial infinity is [H^l 



£0 [1 - 2v^ + (27r - 4qMBH)w^ 



127r 



'(9mbh - 10)^'" + ( — + 4'Zmbh ) + Oiv^) 



(44) 



with 



, 2 -^^mbh'' 



(45) 



3H^orb) 



1/3 



If one transforms to the Jordan 
A'^^/Aq with Aq and a 



and V = (Ml 

frame via TOsco = ^o™sco, r 

given by Eqs. ((28)) and (|29)) . then the above agrees exactly 
with the results of Will and Zaglauer [52| and Ohashi, et 
al. m. 

The fuUy-relativistic, scalar energy flux computed in 
the Teukolsky approach from the numerical solution to 
the equations presented in Sec. IIIII reads 



^ Total 



1=1 m=0 



(46) 



Notice that the I = 1 contribution to the sum docs not 
here vanish, due to the presence of dipole radiation. No- 
tice also that there are two contributions: one that es- 
capes to spatial infinity and one that is absorbed by the 
event horizon of the supermassive black hole. Although 
the latter is usually small relative to the former, and thus 
it is mostly neglected in the literature, we here include 
it. When numerically solving Eq. (|42p for a given Qm-bu, 
we truncate the sum in I when the series evaluated at the 
ISCO, 

''sco = ''iscoj converges to one part in 10^ or bet- 
ter. This requires summing up to Z = 17 for (/mbh = 0.99, 
but only up to Z = 6 for ^mbh = 0. Such a scheme then 
implies that our numerical data is accurate to one part 
in 10'"', which is sufficient for this study. 




S Inf, Wgjj=4e4 
SInf.BD, Wgp=l 
GWInf 



0.0001 



0.001 
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FIG. 2: PN expansion of the (j/^-normalized) gravitational 
wave energy flux (black solid curve) and the scalar energy 
fiux carried out to infinity for Brans-Dicke theory with ljbd = 
4 X 10" (dashed red curve) and uj-bo = 1 (dotted blue curve) 
as a function of velocity for a supermassive black hole with 
spin q = 0.99 and a small compact object with Sgco = 0.188. 
Observe that the scalar flux dominates over the GR one at 
sufficiently small velocities. 

With this at hand, one might wonder how the gravi- 
tational wave flux compares to this scalar flux. Figure [2] 
plots the PN expansion of the (j/^-normalized) gravita- 
tional wave energy flux (black solid curve) and the scalar 
energy flux carried out to infinity with oj^d = 4 x 10^ 
(dashed red curve) and Wbd = 1 (dotted blue curve) 
as a function of velocity for a supermassive black hole 
with spin (/mbh = 0.99. Observe that the scalar flux 
becomes larger than the GR gravitational wave flux at 
sufficiently small velocities. Saturating the Solar system 
bound {oJbd > 4 X 10^), the gravitational wave energy 
flux dominates for velocities v > Wth, where wth ^ 10""^. 
Notice, however, that the difference between the gravita- 
tional wave and scalar fluxes is smallest close to vth- That 
is, gravitational waves in massless scalar-tensor theories 
are the most different from GR waves for v < wth, far 
away from the ISCO. 

Let us now compare the analytic expression for the 
scalar flux to the numerical one. Figure [3] shows the nu- 
merical energy flux - normalized by (av)^ - carried to in- 
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■ inf, Newtonian 
inf, PN 






FIG. 3: Scalar flux at infinity (black solid curve) and at the horizon (blue dashed curve) for a generic massless scalar-tensor 
theory, normalized by (au)^ as functions of the velocity v for Hs = 0, Qmbh ~ 0.99 (left panel), Qmbh = 0.9 (middle panel) and 
Qmbh = 0.6 (right panel). For comparison, we show the PN formula (|44|l (green dot-dashed line) and its truncation at — IPN 
order (red dotted line). 



finity j^^>n='==i<=== (black solid) and absorbed by the super- 
massive black hole jr^'^''^^'''^^ (dashed blue) , as well as its 
leading-order (Newtonian) approximation T^^l'^t"'' (dot- 
ted red) [13 and its higher-order PN expansion T^'^p^"'"'" 
(dot-dashed green) [5Q| as a function of velocity. First, 
observe that the flux absorbed by the horizon is always at 
least an order of magnitude smaller than that carried out 
to infinity. Second, observe that as the small compact 
object approaches the ISCO, the full numerical results 
deviate from the analytic approximations. The disagree- 
ment is larger in the near-extremal case, although even 
then the 2.5PN approximation of Eq. pi)) does remark- 
ably well. For non-spinning or slowly-rotating supermas- 
sive black holes, even the leading order approximation 
reproduces the exact results to high accuracy. 

For corotating orbits, the effect of the high-order PN 
terms absent from Eq. (j44p is to reduce the overall mag- 
nitude of the flux, so that the PN formula overestimates 
the exact flux. In fact, the leading-order, Newtonian flux 
formula can overestimate the flux by more than an order 
of magnitude for orbits close to the ISCO of a maximally 
spinning black hole. Moreover, for highly spinning black 
holes, the horizon flux also contributes to reduce the to- 
tal flux (for Qmbh ^ 0.36, any orbit outside the ISCO 
satisfies the superradiant condition, so that the scalar 
flux at the horizon is negative, cf. Fig. [7] later presented). 
On the other hand, for counter-rotating orbits, the nu- 
merical flux tends to be larger than the corresponding 
PN formula. Since in this case the ISCO is located far 
away from the horizon (6Mmbh < ''isco < SA/mbh for 
< Qmbh < 1), the deviations from the PN formula are 
smaller. 

Our numerical results can be also used to estimate 
some higher-order PN terms missing from Eq. We 
have considered the following corrections to the expres- 
sion in square brackets in Eq. (j44]): 

+03 log (w)] + [04 + as^MBH + a^qliBB 

log(^;)]^;^ (47) 





Fit I 


Fit II 


ao 


12.613 


0.15240 


ai 


-30.520 


-19.619 


02 


13.742 


7.8369 


as 


-1.2497 


-6.5969 


04 





16.569 


as 





-21.100 


ae 





13.959 


a^ 





-1.4791 


ag 








ag 








1-7?^ 


6 X 10"" 


3 X IQ-'* 



TABLE I: Coefficients of the extra term (f47|l in the dipole 
formula (|44|) obtained by fitting our numerical flux. As a 
measure of the fit error, we present the deviation of the coef- 
ficient of determination, E? , from unity. 



The structure of this fitting function is based on the 
3PN and 3.5PN structure of the energy flux expan- 
sion in GR. The coefficients a; can be computed by fit- 
ting our numerical data, obtained spanning the region 
^MBH e [-0.99,0.99] and v G [0.01, Wisco]- We attempted 
approximately a hundred different fits, varying the num- 
ber of fitting parameters, the number of PN order cor- 
rections, and the accuracy of the data to be fitted. The 
best fits we found, ie. those that provide the smallest 
errors with the least number of additional parameters, 
arc shown in Table |T] for two different fits. Notice that 
these parameters arc not meant to represent the next or- 
der terms in the PN scries, but they are just a fit to our 
numerical data to improve the analytical representation 
of the function. 

Figure |4] shows the fractional difference between the 
numerical flux and Fit I (blue dashed curve) or Fit 
II (dotted red curve), which should be compared with 
the difference between the numerical flux and the bare 
PN result of Eq. (|44|) (black sohd curve). This figure 
presents results for a supcrmassive black hole with spin 
(Zmbh = 0.99 as a function of velocity. While the PN for- 
mula can deviate by a factor 300% close to the ISCO, Fits 
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I and II introduces errors smaller than 1% and 0.02%, re- 
spectively. Numerical inaccuracy of our Teukolsky-based 
flux at small velocity {v < 0.1) prevents us from obtain- 
ing reliable fits that deviate less than 1 part in 10'''. Note 
that the fractional difference shown in Fig. |4] represents 
an upper bound; the errors of our fits decrease for less 
extreme values of (/mbh- This is because the PN flux of 
Eq. (j44p agrees better with the numerical flux for lower 
values of ^mbh and because the ISCO velocity decreases 
with decreasing ^mbh- 



le+OO 
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le-06 
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FIG. 4: Fractional difference between the fits, jra^""""!""" ^ with 
parameters from Table |I] and Teukolsky-based fluxes , J-rcuk, 
as a function of velocity for spin 5mbh = 0.99 and compared 
with the fractional difference with the PN flux (|44|) . These 
curves are an upper bound of the flt error: the fractional 
difference for Fit I, Fit II and for the PN formula decreases 
for smaller values of Qmbh. 

Finally, let us consider the effect of this higher-order 
PN terms in the flux on the Fourier transform of the 
response function in the stationary-phase approximation. 
That is, we solve a generalized Fourier integral assuming 
that the gravitational phase oscillates much more rapidly 
than its amplitude (see eg. [iO, HH). After averaging 
over beam-pattern functions, the Fourier transform of 
the response can be written as 



h{f) = A'——r 



7/6gi*(/) 



(48) 



where ^ is a constant amplitude, M = v^l^M is the 
chirp mass, is the luminosity distance, and / is the 
gravitational wave frequency. The Fourier phase can be 
computed from 



*(/) =27r 



//2 p> 



F' 



dF' . 



(49) 



where f'orb is the rate of change of the orbital fre- 
quency, which can be computed by the chain rule: 
i^^^ = Eb{dEb/dF^^.^)-'^. The quantity is the rate 
of change of the binding energy, which from the bal- 
ance law satisfies E^ = — J^rotai, where J^rotai is the total 



energy flux, while dE^/dF can be obtained by differen- 
tiating the binding energy E},, given later in Eq. (|62p . 
The total energy flux J^Totai = -T^gr + ■FtIZ^i"'^'""' , where 
J'^R = {32/5)iy'^v^°[l + 0{v'^)] is the gravitational energy 
flux in GR, while we will model J^To"ar°'°°° = -7^m"pT'°°''- 
We should note here that wc have neglected the contri- 
bution to the rate of change of the binding energy due to 
horizon absorption effects. These will enter at 0{v^) rela- 
tive to the leading-order, Newtonian effect (see e.g. [6^). 
We here only concentrate on high PN order corrections 
induced by the energy flux carried out to spatial infinity, 
such that horizon absorption effects can be systematically 
included in the future if needed. 

Using this and performing the integration, the Fourier 
phase becomes 5'(/) = ^'gr(/) + where the GR 
Fourier phase is 



27r/t, - 




128 ^ 


/3715 


55 


I,) 1.-2/5^2/3 


I 756 ^ 


~ ~9 



Oiu) 



,(50) 



with {tc,4ic) the time and phase of coalescence and u = 
irMf the reduced gravitational wave frequency. Higher- 
order terms in the GR Fourier phase in the PN approxi- 
mation can be found eg. in . The Brans-Dicke correc- 
tion is 



5 Ss 



^,-7/3^2/5 



3584 (wbd + 2) 



7 



^-2/5^2/3 



350 



35 2 

9 ^MBH g 



MBH I ^ 



(51) 



Notice that the leading-order correction agrees identi- 
cally with the results of Q in the large limit, while 
the higher PN terms are new. Notice also that the Fourier 
phase presents similar features to the energy flux. That 
is, the higher PN order terms seem to counteract the 
leading-order, Newtonian term, leading to a smaller num- 
ber of additional gravitational wave cycles induced by 
the Brans-Dicke correction. This then implies that cal- 
culations obtained by using the leading-order, Newtonian 
term only will tend to overestimate the effect of massless 
Brans-Dicke theory. 

One can now map the deformed Fourier phase to the 
ppE framework (2^ . In the latter, one postulates a wave- 
form family of the form h = h^jiSh, where /igr is the GR 
Fourier transform of the response function, while Shp^^ 
is a parametric deformation. The latter can be written 
as (5ft-ppE = (1 + <5^ppe) cxp(z(5\E'ppe), where SA^^^ and 
(J^fppE are amplitude and phase corrections that depend 



on frequency, 
pand (5AppE = 



In its simplest incarnation, one can ex- 



au 



and (5*, 



/3u , i.e. as a parametric 



power series in u, where {a,f3,a,b) are ppE parameters. 
Glcarly, in our case the simplest realization will not do; 
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one must use a more general parameterization, such as 
Eq. (45) in [22|, namely 5A^^^ = and 

K 



= (f>ku 



(fc-7)/3 



(52) 



To match the high-PN order expansion of the phase in 
Eq. (|5ip. we need five terms with the mapping 



5 (sscc 



3584 (wbd + 2) 
'0 , 



2/5 



</>! =0, 



2 



(57r - lOgMBn) ^'^4>Q , 



35 
"9 
84 

5 



^MBH 



350 



(53) 



We see then that the next-to-simplest incarnation of the 
ppE framework easily manages to map to the high-PN 
order Brans-Dicke prediction for the Fourier phase. 



V. MASSIVE SCALAR ENERGY FLUX 

Let us now concentrate on massive scalar-tensor the- 
ories. First, let us consider the scalar field energy flux 
carried out to infinity by a massive scalar field. Solving 
Eq. (fT6|) in the large distance limit, one obtains 



J"" 



3/2 



(54) 

where H^x) is the Heaviside function and we have only 
accounted for the dominant mode. Notice that this flux is 
non- vanishing only for masses that satisfy fi^jb — /.t^ > 0, 
which implies MmbhMs < For typical EMRI veloci- 
ties, V ^ 0.2, the scalar field energy fiux carried out to 
infinity is then non-vanishing only if /i^ < 10~^/Mmbh- 
For larger masses, the scalar mode cannot overcome the 
massive potential barrier and it becomes trapped. These 
results can be generalized to an arbitrary harmonic mode 
m via J7orb < l^s/'m,^ or equivalently Mmbh/^s < mv'^ . 

Figure [5] shows the scalar flux carried to infinity 
and that absorbed by the supermassive black hole with 
^MBH = 0.9, /isMMBH = 0.2 and for different values of I at 
which the series of Eq. P5|) is truncated. Observe that 
the contribution to the scalar flux at inflnity vanishes at 
the orbital distance that corresponds to r2o,-b = l^s/fn, 
in agreement with Eq. ([M)) . Hence, at sufficiently large 
orbital separation, the contribution from the flux at the 
horizon will eventually become dominant, as it occurs 
at w < 0.45 in Fig. [SJ This figure also shows the rate 
of convergence of the scalar energy flux as a function of 
number of I modes summed over. Observe that although 



the flux absorbed by the horizon converges rather fast, 
the one carried out to infinity requires a large number of 
I modes. Overall, however, the total massive scalra flux 
is generically smaller than the massless scalar flux (com- 
pare Fig. [5] with the middle panel of Fig. [3]). If this were 
the whole story, we would therefore generically expect 
smaller scalar-tensor corrections in the massive case. 



le-03 rr 




FIG. 5: Scalar energy flux, normalized by (au)'^, carried to in- 
finity and into the supermassive black hole horizon as a func- 
tion of the velocity for HsM-mbh = 0.2 and Qmbh ~ 0.9. Dif- 
ferent curves are obtained by truncating the series in Eq. (|46|) 
at different values of /. Observe that below a certain velocity, 
^ 0.4, the contribution to infinity vanishes in accordance 
with Eq. (|54p . The resonance corresponding to / = m = 1 is 
also shown. 



Massive scalar-tensor theories, however, have the most 
interesting feature of generically leading to certain res- 
onances in the part of the scalar flux that impinges on 
the supermassive black hole horizon, as recently studied 
in [3^. Figure [S] shows such a resonance at i; '--^ 0.5838. 
These resonances occur when the orbital frequency ex- 
cites the quasinormal frequencies of oscillation of the 
scalar field, which ar e ap proximately proportional to the 
scalar field mass [63l - [65ll. In the limit fj.s <C 1/A/mbh, the 
resonant frequency is [58|: 



1 - 



/^S-^^MBH 

l + l + n 



-,1/2 



0,1, 



(55) 



Clearly then, rin^g ~ /is in the /Xs-M^mbh <C 1 limit. We 
can relate this to the orbital radius via 



yMBH 
(A/MBHO„b)2/3 



(56) 



where we define the resonance location as rR^g = r(D,j^^^). 
Obviously, these resonances are not present in the mass- 
less case as then n^^, — >• 0. 

Such resonances greatly enhance the magnitude of 
■Z^r™'°° and, in turn, the amount of energy flux radi- 
ated [39|. In Appendix El we compute the resonant flux 
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analytically in the fisM^BH ^ 1 limit, for any {l,m) 
and for the fundamental mode (cf. Eq. (|A17p ). When 
I = rn = 1, Eq. ()A17P reduces to 
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2 -^^MBH 
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(57) 



where B = l+AK.^ and /C = -2A/MBHr+fc+/(r+-r_). We 
have found very good agreement between this expression 
and our numerical solutions for /ZsMmbh ^ O-l- 

In the small /is limit, the resonance location scales as 
TRes ~ Ms ^^"^ ^iid the scalar flux at resonance grows in 
magnitude when fig 0. For generic / and Qmbh 0, 
the peak flux scales as -7^+ peak™ ^ [cf. Eq. (|A17|) 

in Appendix ((^ ]. For very small ^mbh, the peak flux 
at resonance is instead positive, and it can also be very 
large. For example for a Schwarzschild black hole, 
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(58) 



In the small /ig limit, the resonance width is twice the 
imaginary part of the quasinormal mode frequency [58| 
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(59) 



For the fundamental mode, n = with I = m = 1, and 
in the small /is limit, the equation above reduces to 

1 



12M„ 



(/isMMBn) (^MBH - 2r+/is) . (60) 



Notice that MMBnAfi cx (/IsMmbh)^ and it is thus in- 
credibly small for small /igMMBH- For larger values of 
/<sMmbh, the width can be computed numerically, using 
a continued fraction method (63l - [65j . Some numerical re- 
sults are shown in Table HIl where Ai? = Afldr/dil is 
the resonance width evaluated at r = Vj^^^ using Eq. (|56p . 



Scalar flux resonances can lead to certain floating or- 
bits, as recently studied in [s^. Such floating orbits are 
defined as those where the scalar flux identically can- 
cels the 0{i'^) gravitational wave flux, thus greatly slow- 
ing down the inspiral. The inspiral evolution is then 
sourced by the shedding of mass and angular momentum 
of the supermassive black hole. Cardoso et al. [Ill found 
that such orbits are possible outside the ISCO, provided 
qmbh > 0.36 and the scalar mass /is is small enough. Sec- 
tion IVIII will consider how such resonances can be used 
to constrain massive scalar-tensor theories. 
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AR/(2Mmbh) 


0.35 


1.52818949793075 


-0.0522 


1.15 X 10" 


-V 


0.3 


1.78503938021340 


-0.0752 


1.52 X 10" 


-7 


0.25 


2.10016240393323 


-0.0969 


5.93 X 10" 


-8 


0.2 


2.53500275855866 


-0.1200 


1.82 X 10" 


-8 


0.15 


3.18939434130550 


-0.1467 


3.30 X 10" 


-9 


0.1 


4.33400288873563 


-0.1828 


4.38 X 10" 
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TABLE II: Orbital radius at resonance and peak scalar flux 
for n = 0, Z = m = 1, Qmbh = 0.99 and several values of 
fisMijiBH- Ai? is the resonance width in the radial dimension. 



VI. PROJECTED CONSTRAINTS ON 
MASSLESS SCALAR-TENSOR THEORIES 

In this section we will study the constraints one could 
place on massless scalar-tensor theories given a gravi- 
tational wave observation of an EMRI. We begin with 
a review of the EOB approach [H, 113, [66I - I8H . fol- 
lowing in p articular the model developed by Yunes, et 
al. [H, l36l [86| . We then continue with a discussion of 
how this model can be modifled to include scalar-tensor 
modiflcations. We conclude with a numerical evolution of 
orbits in the EOB-EMRI scheme to determine plausible 
projected constraints. Note that, when needed, we shall 
map the metric perturbations from the Einstein frame to 
the Jordan frame, by using Eq. (j23p . 



A. EOB Modeling of EMRIs in General Relativity 

The use of the EOB-EMRI approximation is not funda- 
mental to this paper. First, the only EMRI orbital regime 
of interest for data analysis is the inspiral, as almost all of 
the signal-to-noise ratio builds during this regime. There- 
fore, the EOB-EMRI approach does not account for the 
plunge and merger, as these phases are irrelevant for EM- 
RIs. Second, we choose an EOB-EMRI model because 
modifying it to include scalar-tensor theory corrections 
is fast and straightforward. We could have instead used 
a time-domain Teukolsky evolution [13, [13 , or an inter- 
polated frequency-domain Teukolsky evolution [H, [o^l ■ 
However, in Refs. [H,[H,|1^, it has been shown that an 
EOB-EMRI approach is equally valid during the inspi- 
ral phase up to the ISCO for quasi-circular orbits in the 
equatorial plane. 

The EOB-EMRI framework makes use of the so-called 
adiabatic approximation, i.e. the radiation-reaction time- 
scale is much longer than the orbital one. In this approx- 
imation, the gravitational wave phase can be obtained by 
solving 



Ei,{n„^), </) = ri„., (61) 



where is obviously the small compact object's orbital 
angular frequency and (j) is the orbital phase, is the 
system's binding energy and Eb = — J^rotai is the total 
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rate of change of this energy, where J"Totai is the energy 
flux. 

We model the binding energy of the system as that of 
a test-particle in orbit at radius r around a Kerr black 
hole [ill: 



Eh 



1 2Af"]viBH/^SCO i ^MBH -^^MBh/^SCO 
1 3A'/mBh/^SCO i 2^]viBH -^'-^MBh/''"sCO 



(62) 

where the ± stands for prograde or retrograde orbits. 
This is the exact Hamiltonian for a test-particle in orbit 
around a Kerr black hole (without the constant rest mass 
piece). Therefore, Hamilton's equations are equivalent to 
the evolution of geodesies in the Kerr spacetime. 

In GR, the only energy sink is that induced by the 
emission of gravitational waves. Because of this, one as- 
sumes the balance law Ejj = — J^cR.Totaii where the lat- 
ter is the total energy flux J"GR,Totai = -^gr,oo + -^gr,+ - 
The component that describes radiation that escapes 
to spatial infinity is labeled J^gr,ooj while the one that 
describes radiation absorbed by the supermassive black 
hole is labeled T^-^ ^. The former can be approximated 
as [Zl [tI HI 



8 i 



(63) 



1=2 m=0 



where the multipole-decomposed gravitational wave can 
be re-expressed in the following factorized form 



1 NGwt,ep Qtp rp 



{PlniY 



(64) 



Here, tp is the parity of the waveform {i.e., = if 
£ -f m is even, ep = 1 if £ -I- to is odd), while Sy (v\ 
Tim{v), 5im{v) and pem{v) can be found in [lllzl,!!!!. 
The Newtonian waveform is simply 



MBH (Cp) 



R 



where 



(65) 

are spherical harmonic functions, while 



and C£+ep are numerical coefficients [73 



The fiux presented above only accounts for the amount 
of radiation that escapes to infinity -^gr.oo, to which one 
must add the amount of radiation absorbed by the super- 
massive black hole J-^^ ^ . In the EOB-EMRI framework, 
this is accounted for by linearly adding black hole absorp- 
tion terms and calibration coefficients, fitted to a more 
accurate, numerical fiux [36j . The calibration accounts 
for unknown, high-order PN terms. 

With all of this at hand, one can solve the differen- 
tial system in Eq. (|61[) with certain post-circular initial 
conditions [131 • We implement this by first perform- 
ing a mock-evolution that starts at IOOMmbh with post- 
circular initial conditions (see e.g. [s^) and then reading 
off [0(0), 0(0)] at the desired starting point of the real 
evolution. The gravitational wave phase and amplitude 



can be obtained from the waveforms via 

hpr, 
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(66) 



Equations ((62)) and (|63)) neglect the effect of the small 
compact object on its own evolution, i.e. the correction 
to the metric, Hamiltonian and radiation-reaction due to 
the small compact object's self-gravity. These correspond 
to conservative and second-order dissipative self-force ef- 
fects, which we will here neglect. Although such effects 
are important for the modeling of EMRI waveforms, they 
are not currently known in the strong-field regime and 
thus cannot be properly incorporated into this framework 
(however, see the recent work [92j where such effects have 
been included in the case of Schwarzschild black holes 
within GR). 



B. Scalar-Tensor Modification 
to the EOB-EMRI Model 



The EOB-EMRI framework described in the previous 
subsection can be modified to model EMRIs in scalar- 
tensor theories by properly accounting for all energy 
sinks. In scalar-tensor theories, the gravitational wave 
energy flux is complemented by the scalar flcld energy 
flux: 



^ GR.OO 



GR,+ 



(67) 



Depending on the particular scalar-tensor theory stud- 
ied (eg. massive or massless), may or may not con- 
tain resonances. Massless (^s = 0) scalar-tensor theories 
are considered in Sec. IVI C| by modeling the energy flux 
via the numerical results presented in Sec. I IVI The de- 
tectability of massive scalar-tensor theories is discussed 
in Sec. IVIII Of course, the modification of the energy flux 
due to scalar emission will induce a non-GR radiation- 
reaction force that, although small at any particular time 
(a part from possible resonant effects), will build up pos- 
sibly leading to strong modiflcations in the waveforms. 

In addition to this, there is in principle a modification 
to the binding energy. This effect, however, is difficult 
or impossible to measure with gravitational waves. This 
is because in Brans-Dicke theory, the scalar field can be 
thought of as renormalizing Newton's constant, or equiv- 
alently the bare system parameters. For example, given a 
gravitational wave detection, one would think one is mea- 
suring the bare component masses, while in reality one 
would be measuring a certain component tensor mass, 
composed of the product of the inertial masses and the 
Brans-Dicke coupling parameter [l^, . For this rea- 
sons, we will ignore these effects in the current paper, 
and only model dipolar corrections due to scalar field 



emission. 
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C. Projected Constraints on Massless 
Scalar- Tensor Theories 

Consider the following two systems: 

• System F: A/mbh = 10^ Mq, TOsco = 1.4 Mq (i/ = 
1.4 X 10-7), (Zmbh = 0.99, t/Mmbh e (1.927, 1.455), 

V € (0.649,0.714), few & (0.00176,0.00235) Hz, 
total cycles in ^ 3.75 x 10^ rads of gravitational 
wave phase. 

• System R: Mmbh = 10^ Mq, ttisco = 1-4 Af© (i/ 
1.4x 10-5), g„3H = 0.7, r/MMBH G (18.145, 15.985), 

V e (0.234,0.249), few e (0.008,0.010) Hz, total 
cycles in ^ 1.79 x 10^ rads of gravitational wave 
phase. 

The one-year evolution of System F within a LISA-like 
band samples a strong-field region and it is representa- 
tive for a configuration which maximizes the relativistic 
effects of the scalar emission. The small compact ob- 
ject orbits very close to the supermassive black hole and 
indeed the evolution is stopped at the orbital separation 
corresponding to the ISCO. On the other hand. System R 
is a more realistic configuration, in which the supermas- 
sive black hole is rapidly spinning, but is not extremal, 
while the orbital velocity and the mass ratio are optimal 
for the new eLISA configuration [4l|. In this case, the 
EMRI samples a moderately strong-field regime and the 
evolution is stopped at the orbital radius corresponding 
to a final gravitational wave frequency of O.OlHz. 

The calculation of the gravitational wave phase in 
massless scalar-tensor theories will depend not only on 
the value of a, or equivalently Wbd, but also on the sen- 
sitivity parameter Sgco. As for the former, we will here 
work directly in terms of Wbd and we will saturate the 
Solar System constraint ojbd = 4 x 10*. This will allows 
us to determine whether such a value of cJbd leads to an 
observable effect in EMRI waveforms. As for Sscoi we use 
the representative value Ssco ~ 0.188, which corresponds 
to the average of three sensitivities for three different neu- 
tron stars of the same mass but with different equations 
of state [5^ . 

Before comparing waveforms directly, we will perform 
certain shifts in time and phase so that the waveforms 
are initially aligned in phase and frequency. In partic- 
ular, we minimize the statistic in Eq. (23) of Ref. [75| . 
which is equivalent to maximizing the fitting factor over 
time and phase of coalescence in a matched filtering cal- 
culation with white noise [t^. We carry this out in the 
low-frequency regime, inside the time interval (0, 64)Agw 
(roughly 0.01 months), where Aqw is the gravitational 
wave wavelength). 

Let us now compare the phase and the amplitude of 
GR and massless scalar-tensor theory waveforms. We 
define the dephasing by Scj) = (/)gr — 0st, measured in 
radians, where 0gr is the waveform phase in GR and 
(/)sT is the waveform phase in scalar-tensor theories. Wc 
define the normalized amplitude difference by 5 A = |1 — 



Ast/^gr|, where Agr is the GR waveform amplitude and 
is the waveform amplitude in scalar-tensor theories. 
Figure [S] shows the dephasing and normalized amplitude 
difference as a function of time in units of months for 
Systems F (left) and R (right), respectively. We show 
results obtained using both our numerical flux and the 
FN formula in Eq. pi)) . 

Our results show that EMRIs cannot place constrains 
on massless scalar-tensor theories, including Brans-Dicke 
theory, that arc stronger than current Solar System 
ones. Figure [S] shows that, even considering one of the 
most relativistic cases among detectable EMRI systems, 
the accumulated dephasing over one year is less than 
Scf) ^ 0.1, which is roughly expected to be the limit 
for detection with eLISA (unless the signal-to-noise ra- 
tio is much larger than ten). Of course, if we had cho- 
sen Wbd < 4 X 10'^, then the dephasing would have been 
larger, as the emitted flux scales as but such val- 

ues of Wbd arc already ruled out by Solar System ex- 
periments. Furthermore, the FN results obtained from 
Eq. pi)) overestimate the dephasing and the normalized 
amplitude difference computed with our full numerical 
flux. This implies that tests of Brans-Dicke theory with 
EMRIs will be more difficult than previously expected 
by extrapolating leading-order, FN results to the strong- 
field region. 

The right-panel panel of Fig. [5] shows that our results 
perfectly agree with the FN expectation when consider- 
ing System R. Indeed, System R is less relativistic than 
System F, and thus, for the former the 2.5FN formula 
of Eq. is a very good approximation (cf. Fig [3]). 

Furthermore, in this case the dephasing after one year 
evolution is larger than that for System F. This is not 
only because System R has accumulated ten times more 
phase than System F in the same observation time (be- 
cause the radiation-reaction time scale is longer), but also 
because Brans-Dicke theory leads to larger relative devi- 
ations when the EMRI is farther away from the ISCO. 
As we showed in Fig. [51 the massless scalar and GR en- 
ergy fiuxes are of the same order of magnitude for orbital 
velocities v ~ Wth, and thus, the radiation- reaction force 
is most modified for System R as z;th < ^Sys r < "^Sys f. 
This suggests that to get the strongest constraints, one 
should consider EMRIs with v ^ vth, which translates 
to separations around rth ^ vf^^ ^ 10^ A/mbh . Unfortu- 
nately, the larger the orbital separation, the lower the 
gravitational wave frequency, and for such large separa- 
tions, the EMRI emits waves completely outside of the 
sensitivity band of future detectors. 

These results go against the expectation that EM- 
RIs will be the best probes of modified gravity the- 
ories. This expectation derives from the observation 
that EMRIs sample the near-horizon region, and thus, 
the strong-curvature regime of the underlying gravity 
theory [2J|. However, although Brans-Dicke gravity is 
among the most popular alternatives to Einstein's the- 
ory, it is not, by itself, a strong-curvature modification of 
GR. As Fig. [2] showed, Brans-Dicke theory leads correc- 
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FIG. 6: Dephasing and normalized amplitude difference as a function of time in units of months for System F (left) and 
System R (right). The black and red curves are obtained by using the numerical fluxes constructed in Sect. IIIII and their 
PN approximation in Eq. (|44[) . respectively. We consider a Brans-Dicke parameter which saturates the Solar System bound, 
cjbd = 4 X 10*. 



tions that dominate over GR at low inspiral velocities, 
not close to the horizon. EMRIs should experience much 
larger deviations in theories that introduce true, strong- 
curvature corrections, like e.g. Chern-Simons gravity [93| 
or alternative theories with generic quadratic curvature 
corrections [l^]- These theories are poorly constrained 
in the weak-field regime [l^, [1^, while EMRI gravita- 
tional waves are expected to lead to the strongest con- 
straints [lT|.[l3H96tf97t. 



VII. PROJECTED CONSTRAINS ON MASSIVE 
SCALAR-TENSOR THEORIES 



Let us now consider the effect of massive scalar-tensor 
theories on the gravitational wave phase, and try to es- 
timate whether this is observable with future detectors. 
As we noted in Sec. |Vl the total flux in the massive case 
is generically smaller than the one in the massless case. 
Therefore, except for resonances in the scalar flux [39} . 
the dephasing induced in the massive case is smaller than 
in the massless case. Resonances, on the other hand, can 
lead to large dephasings, as we will show in this section. 

We begin by considering the region of phase space 
where resonances are present. We then summarize an an- 
alytical prescription to calculate the gravitational wave 
phase in GR. We choose such an analytical route, instead 
of an EOB-EMRI one, because the resonances tend to be 
very difficult to resolve numerically. With this in hand, 
we then proceed to calculate the accumulated phase when 
an EMRI traverses a resonance and the corresponding 
dephasing. 



A. Floating versus Non-Floating Orbits 

Depending on the spin of the massive black hole, the 
resonances in the scalar horizon flux can lead to two 
qualitatively different effects. As shown in Fig. [71 if 
Qmbh > 0.36 any orbital frequency < fiisco('7MBH) 
also satisfies the superradiant condition for m = 1, 
mbh)- Thus, possible resonant fiuxes are neg- 
ative in this region. On the other hand, if Qmbh ^ 0.36, 
non-superradiant frequencies are also excited and the res- 
onant fiux becomes positive. 

A negative flux leads to the interesting possibility of 
a floating orbit, i.e. an orbit where J^gr + J^'' = 
[331. Recall, however, that this does not imply the orbit 
freezes at a given radius, but rather that its inspiral slows 
down, as it is not driven by the shedding of mass and spin 
angular momentum of the supermassive black hole. The 
floating orbit condition on the fluxes defines the region 
in {^s, ot) space where floating can occur: 

a > "c = ^^^V^^MBH - 2?/(l + 9mbh)2/"^^ ■ 

V 15 

(68) 

where we have assumed y ~ /^s7\/mbh ^ 1 and recall 
that, by Eq. (|15p . the leading-order quadrupolar flux 
reads J^|r = (32/5):/^z;^°. Equation (|68p ensures that 
fc+ < at the floating frequency ri„,b = refloat- No- 
tice that the floating frequency is smaller, but very close 
to, the resonance frequency, with differences of 0{^^). 
This, however, does not imply that T^"""^"""' = 
J^^'""''""™ (refloat), since jg ^ very steep function 

of il close to resonance. 

In Fig. [SI we plot ac in Eq. as a function of the 
scalar mass, /is, in units of 10~^® eV for two systems 
where the supermassive black hole has mass IO^Mq and 
IO^Mq and different spins. Notice that ac is essentially 
insensitive to the spin of the supermassive black hole. 
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FIG. 7: ISCO frequency and angular velocity at the horizon 
as functions of the spin parameter for a Kerr black hole. In 
massive scalar tensor theories, resonances occur for particles 
in circular orbits at florb ~ /is + 0{jj.s). Stability of the orbit 
implies fiorb < i^isco, i.e. the relevant region is below the 
blue line. For Qmbh 0.36, these frequencies always satisfy 
the superradiant condition for m = 1, Qarb < For gMBH < 
0.36, the resonant frequency may also be non-superradiant in 
the region limited by the blue line and the red dashed line. 
Note that resonances may occur even for Schwarzschild black 
holes. 
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FIG. 8: Critical value, Oc, as a function of the scalar mass 
fis, for several values of Qmbh and for Mmbh = W^Mq (lower 
curves) and Mmbh = lO' Mq (upper curves). Floating orbit 
exists for those values of a above the curves and for q > 0.36. 
Curves terminate when Hs ~ fiisco. 



the definition in Eq. (1291) of lObo in terms of a, we find 

3 15 / iV 

Wbd < t^BD = ^2 ^ 16 V ^ ^/ ^ 

^ ^ , 

(Jmbh - 2(1 + Vl - q'i.j,^) Ids Mmbh 

(69) 

The regions in (uj^m fJ-s) phase space where floating oc- 
curs is shown in Fig. [T] for a supermassive black hole 
with A^MBH = lO'''i\/0, spin (/mbh = 0.9 and a neutron 
star sensitivity of Ssco = 0.188. The red-dashed curve is 
the boundary between floating and non-floating orbits for 
this system, given by Eq. ([69]): non-floating resonances 
exist for values of (wbdjMs) above this fine, while float- 
ing resonances exist for values below it. For comparison, 
we plot the (wbd,Ms) region (below the solid black line) 
that is ruled out by the tracking of the Cassini space- 
craft [13, HI]. Notice that as fXs 0, one recovers the 
massless Brans-Dicke bound Wbd > 4 x 10*. Since the res- 
onance frequency scales as rinos ^ Ms (cf- Eq. ([M)) ). we 
also show the values of fig (green, dot-dashed curve) for 
which the resonant frequency occurs inside the classic- 
LISA sensitivity band [(10-^,10-2) Hz]. Finally, the 
vertical dotted line denotes the ISCO frequency, riisco- 
Curves in Fig. ([T|) terminate approximately when ^ 
f^isco, which corresponds to the small compact object 
reaching the ISCO. 

Just because a floating or non-floating resonance ex- 
ists does not necessarily mean that it will lead to a de- 
tectable effect in the waveform. Whether such resonances 
will matter or not depends both on the height of the 
resonance and the amount of time it is active. In the 
non-floating case, this lifetime is just the width of the 
resonance Str- In the floating case, the resonance width 
does not come into play, as floating orbits are slowed 
down and inspiral due to the supermassive black hole's 
backreaction. In the latter, one must wait a time St until 
the supermassive black hole loses enough mass and spin 
angular momentum so that the height of the floating res- 
onance is decreased, allowing the small object to pass 
through. In what follows we will study how both types 
of resonances affect the waveform and whether they lead 
to observable effects. 



B. Analytical Modeling of the Gravitational Wave 
Phase in General Relativity 



while for different symmetric mass ratios i^, ac shifts hor- 
izontally. 

In the particular case of massive Brans-Dicke theory, 
we can also compute the region of the (wbdiMs) param- 
eter space where floating can occur. Using Eq. (p5|) and 



Let us first consider two of the most popular methods 
to compute the gravitational wave phase in GR. The first 
method consists of employing the balance law 

Er = -nn, (70) 

where Ef, is the time derivative of the binding energy and 
J-Qn is the total fiux of energy carried away from the bi- 
nary system. The binding energy can be computed from 
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the Hamiltonian of the system, and to leading, Newto- 
nian order it is simply 



-7: 



(71) 



where we have defined 7 = Mmbh/''- Through the Hamil- 
tonian, one can also establish Kepler's third law, namely 
that 



7 = X = (i\fMBHf^orb) 



2/3 



(72) 



which then leads to E^^- ~ — i(i^A/MBH)/2. The total en- 
ergy flux in GR is simply given by the total gravitational 
wave energy flux, which to leading-Newtonian order is 
simply 



GR 5 



(73) 



where again we have used Kepler's law to relate 7 to x. 

With this at hand, we can use the balance equa- 
tion ([70| to derive a differential equation for the fre- 
quency evolution, namely 



5 A/mbh 

or in terms of the orbital frequency fiorb = 27rfo,b, 



(2^MMBHi^o,b) 



11/3 



(75) 



These differential equations are separable, and thus, we 
can solve them to obtain 



x^-{Q) = \q-^'\ 
4-3/2 

ZTTMmbh 

where we have defined the quantity 



5AfM BH 



(tc - t) 



(76) 
(77) 

(78) 



The quantity tc is a constant of integration chosen such 
that t = tc when F^^^, = +00. This constant is sometimes 
referred to as the time of coalescence, and thus, tc^t > 0. 

With such a frequency evolution, it is simple to obtain 
the GW phase. For the dominant harmonic mode in a 
quasi-circular inspiral, this phase is related to the orbital 
one via = 20„t,, and thus 



47r / F„,(t') dt' 



10 



a;3/2(e') de'. (79) 



This expression agrees identically, to leading, Newtonian 
order, with PN results [sl]. Inserting typical values, we 
find 



GR 
GW 



(tc - t) 



= 4 X lO^'rads 
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(81) 



The GW phase in Eq. (|M)) is exact to leading, Newto- 
nian order, but it might not be very useful for EMRIs, 
since it assumes that the orbiting body reaches merger in 
a one year evolution. EMRIs, however, can be in orbit for 
tens of years if their mass ratio is small enough. When 
dealing with EMRIs that might not merge in a one- year 
evolution, it is perhaps more appropriate to rederive the 
GW phase in a slightly different way. Consider then the 
following expression for the GW phase 



'^GW 



dr 



(82) 



where (?";„;, rfi„) are the initial and final orbital radius 

(74) 

that are observed, while 57„b — 7'^^^/AfMBH by Kepler's 
third law. with the rate of radial inspiral obtained 
by the chain rule: 



dE'i;'^ 



dr 



dt 
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The GW phase in GR is then simply 



,GR 
"^GW 



1 

16^ 



5/2 



5/2' 



(83) 



(84) 



The GW expression in Eq. (1841) can be converted into a 
function of time by solving the radial evolution equation. 
The solution to Eq. ([55]) is simply 



ry-iGR 
-'tin.ini 



where Tfi, 

Tini to Tfi,, 
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256^ 



( ''ini 



(85) 



i„i is the time for the small object to go from 
Using this in Eq. (|84)) . we obtain 



iGR 
^GW 
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(86) 



We can solve this integral to find 



x-5/2 9 

•PgwI^J 161/ V ■ 



(80) 



The time T'f^'^^j can be required to equal the observation 
time, i.e. 1 year. This expression agrees with that used, 
e.g. in [99l. llOOt. and it reduces to Eq. dUl) when the 
second term in square brackets dominates. 
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C. Effect of Non-Floating Resonances 
on the Gravitational Wave Phase 

Let us now concentrate on the gravitational wave phase 
as an EMRI crosses a non-floating resonances. We model 
the correction to the GR gravitational wave flux via T = 
7"^^ + J"', with 



s, massive 



(87) 



where is the height of the resonance, t^^^ is the time 
at which the resonance occurs, 61^^^^ is its width, where 
recall that l-i{x) is the Heaviside function. Notice that 
Eq. ([87)1 is nothing but a top- hat function of width (S<Roa, 
centered at f = t^^^ and with a time-dependent height 
h{€). The height h{i) represents the total scalar flux 
•^To™'"™! which must satisfy < |J-"GR(t)| for non- 

floating resonances. If /i(t) > 0, one has a sinking res- 
onance, i.e. a resonance whose modification to the fiux 
enhances the inspiral rate. Of course, one can also have 
h{t) < but still \h{t)\ < IJ'^rI such that the orbit is 
non-floating and non-sinking. 

For massive scalar-tensor theories, we can compute 
i5tR„a explicitly. This width can be obtained from 



\An\ 



(88) 



can treat h{t) = /i as a constant, i.e. h = 

Eq. dST]) with a < Uc- 

Let us now use the methods in Sec. IVIIBI to derive the 
gravitational wave phase as an EMRI traverses a non- 
fioating resonance. Let us then use the balance equation 
to derive a differential equation for the orbital frequency 
evolution: 



Eh -^GR 1 



(92) 



where Tan is the GR gravitational wave energy flux and 
J^" is the scalar energy flux of Eq. ([57]) . Since the Hamil- 
tonian of the system is not modified by the scalar res- 
onance, we still have that Ei, is given by Eq. ([7T]) to 
leading-order and Kepler's third law still holds. 

The evolution equation for the orbital frequency is then 
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(27rA/MBHF„0 



11/3 



27r vMl 



^ (27rA/MBH^^„,b)'/'-F^ (93) 



This is a differential equation for F^^^^ that is not sepa- 
rable due to the term. However, since we are dealing 
with non-floating resonances, we can search for solutions 
-F^orb = -F^o^b + SFarb, whcrc 5F„^ is assumed to be smaller 
than F^^^. Inserting this ansatz into Eq. (|93|) we obtain 
the zeroth-order differential equation 



where Afi is the resonance width in frequency space 
[cf. Eq. ((59|)]. We can model f2„fb at resonance as 



dEb f dEb 

"dT \~dn 



W GR ~ /rcs 
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(dEb\ 
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(89) 



where in the last line we have neglected the scalar flux, 
assuming that in general this is smaller than the gravi- 
tational one. This then leads to 



6t„ 



5 Mme 
96 1/ 



(A/mbh^^RoJ""^' (MMBnAf^) , (90) 



where we have used that J^Ir = {32/5)i'^x^. Substituting 
in An = Ar2'™" from Eq. ([60]) and ^2r„, from Eq. (|55l) . 
we find 



5 A/me 

1152^7 



(^sA/mbh)^'^^^ (^mbh - ^r+fls) , 
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Qmbh 
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(91) 



where in the second equality we have scaled the width 
by typical numbers. Clearly then, St^i^, <C T^i,^ and we 



5Tr ML 



(27rA/MBHF':b) 



■GRai/3 



(94) 



and the first-order differential equation 
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(95) 



The solution to these differential equations is sim- 
ple. The solution to the zeroth-order equation is that 
of Eq. ((77)) . The solution for the correction to the fre- 
quency evolution can be written as 



SF„ 



271 vMl 



^—e^ f (2^A/MBH^^r )'^' e-^T'dt' , 

^MBH J 



(96) 

where we have set the integration constant to zero and 
we have defined the quantity 



X 
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n r-^'^ dx^ 

~2 J ~ 
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(97) 



With this at hand, Eq. (j96p becomes 
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27r i^M^^^ 



x-^{t') ^,{t') dt', 



h (5<Re. , 



(98) 
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where in the second Une we have performed the integral 
over the top-hat function Tg and evaluated all quantities 
ait = , as this is the only time the resonance is active. 

From this expression, we can now find the waveform 
phase. Using Eq. ([7^. the gravitational wave phase is 



67W 



-2, 1/3 U Xf T 

"-Res fJt'Rcs -'-obs : 



(99) 



where we have used that u{f) = nMf, with / the grav- 
itational wave frequency, and the second term integrates 
simply to Tobs because it is constant. The gravitational 
wave dephasing is then simply 



(100) 



Of course, recall that this expression is only valid when 
J^* is assumed much smaller in magnitude than J'q^^. 

The above results are generically applicable to any 
non-floating resonance that is active for a very short 
time, such that their modifications to the flux can be 
treated as effectively a constant small relative to the GR 
flux. In particular, the above dephasing is also applicable 
to th e GR reso nances discovered by Flanagan and Hin- 
derer lOll Il02{ . Requiring that the resonances lead to an 
observable effect (i.e. a dephasing larger than 1 radian), 
the product of the timescale and the magnitude of the 
flux relative to the GR flux must satisfy 
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GR 



^1 u^,', 
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(101) 



where we have rescaled all quantities by typical values. 
The incredible smallness of this dephasing can be traced 
back to the fact that the resonance is extremely narrow 
[cf. Eq. (HI])]. 

We can now require that this dephasing be greater than 
one radian and solve for the value of a that would allow 
for this. In order to have a dephasing larger than one 
radian, a > anf, where the latter is defined as 



anf « 482.9 




0.01 



-8/3 



(104) 



Alternatively, we can use Eq. ([29]) and write the corre- 
sponding bound for Wbd < '^bdi where 



.nf 



10" 



IO^Mq 



0.188 J V lyr 

n,r -\ 16/3 
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(105) 



Such a critical value of is shown in Fig. [T] as a dot- 
dashed, blue curve. That is, for values of (wbdjMs) be- 
low this curve, the dephasing is larger than one radian, 
while for values above the curve the dephasing is less 
than one radian. Observe that the red-dashed curve is 
always above the blue dot-dashed curve, i.e. there are 
no values of (wbdjA's) that can simultaneously lead to a 
significant dephasing (larger than 1 rad) and to a non- 
floating resonance inside the classic-LISA frequency band 
(green, dot-dashed curve). Mathematically, there are no 
values of (a^fis) that can simultaneously satisfy a > a^f 
and a < ac- This then automatically implies that non- 
floating resonances cannot be observed with a LISA-like 
mission. 



D. Effect of Floating Resonances 
on the Gravitational Wave Phase 



Equation (jlOip is an interesting constraint: it implies 
that if the relative height is of 0(10""'^), then the reso- 
nance width must be at least 4.2 seconds wide for the 
corrections to the GR flux to have an observable effect 
on EMRI waveforms. 

Let us now evaluate this dephasing for the height h and 
width (5f:Res appropriate to massive scalar tensor theories. 
Putting all pieces together, the dephasing becomes 
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Evaluating this equation to leading-order in HsMmbh 
1, we obtain 
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Consider now the gravitational wave phase as an EMRI 
traverses a floating resonance. During such resonances 
the total energy flux can be written as 



GR 



(106) 



GR — 0{v^) is the GR gravita- 
tional wave energy flux and = 0{v'^) is the scalar en- 
ergy flux. As we shall see next, if these two terms cancel 
each other, for example at a floating resonance, an ad- 
ditional non-adiabatic contribution, J^non-adiab = 0{iy^), 
arises due to mass and spin angular momentum loss of the 
supermassive black hole. Therefore, at resonance the bi- 
nary system continues to lose energy and inspiral, albeit 
at a slower rate, as the supermassive black hole back- 
reacts by shedding mass and angular momentum. 

We can make this clearer by considering the balance 
law Eh ~ — J^Totai- The binding energy of the small com- 
pact object, Eq. (]H^ . can be written in terms of the 
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orbital frequency as 



1 - 2w2(l - V^Qmbu 



,1/3 



(107) 

where we recall that v = (MMBH^^orb)^^'^- Equation (|107p 
asymptotes to Eb ^ i^M^bh — {'^Mmbh/'2)v^ in the w <C 1 
limit, which agrees with Eq. (|71l) . To leading order in 
I/, Mmbh and (7mbh arc constant, but to next order they 
evolve: 



Eh 



dEi, 



dEb 

dq^Bi 
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where a dot denotes a total derivative. We can then use 
the balance law to solve for namely 

dEb\ ■ ^ f dEb 



0„ 



dAL, 



dEb 



(109) 



which is exact to all orders in v. While the first term on 
the right-hand side is oiO{v^), the second term (in paren- 
thesis) is of 0{v^\ because (-Mmbhj9mbh ) = 0{v^) and 
{dEb/dMMBH,dEb/dqMBH) = 0{y). Usually, one thinks 
of (A/mbh, ^mbh) as constant, and through the chain rule 
identifies the left-hand side as simply Eb. In that case, 
one can effectively think of the term in parenthesis on the 
right-hand side as the additional J"„o„-adiab in Eq. ()106p . 
We emphasize that we are not adding an J^„o„-adiab con- 
tribution to J-Votai, but rather the time-dependence of 
(Mmbh, ^mbh) serves as an effective modification to J^xotai- 
The evolution of the orbital frequency at resonance can 
be computed explicitly by using the appropriate rates of 
change of mass and spin angular momentum. First, we 
recognize that (Jmbh = JmBul M^mBin due to the laws of 
black hole mechanics. Moreover, for circular orbits we 
also have that (JA/mbh = f^orbi^AiBH, and thus 

'^"""^ M — n M — ■ ^^^"^ 

Then, substituting 

M^B^FZT) = K^FZT) « ~n.{FZT) , (111) 



and Eq. pTU)) into Eb = Q with Eq. PU5)) . and solving 
for f^orb we obtain [33, 



r2o,b ^ 32i/^ M, 



2 A/A^/^. rj^^/^ 



(112) 



in the small /i^ limit. Therefore, once an EMRI hits a 
floating radius, its inspiral evolution slows down, as de- 
scribed by Eq. (|112p . which is 0{v) smaller than the GR 
evolution of the orbital frequency. This evolution equa- 
tion kicks in after the particle has hit a floating resonance 
and it lasts for a time 5t. 



The floating resonance lifetime 5t is defined as the time 
that it takes the supermassive black hole to lose enough 
mass and angular momentum for the scalar flux to not 
be able to compensate the gravitational flux, and thus, 
5t is controlled by J^,on-adiab- The floating lifetime can be 
computed from 



5t = E- 
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-|-,pcak 
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The rate of change in the resonant peak can be com- 
puted by differentiating Eq. ([57]) (considered as a func- 
tion of riorb, Mmbh and (7mbh), applying the chain rule 
as in Eq. (|108p and then using the rates of change Ct^^t, 
Mmbh and (7mbh computed above. Finally, in the small 
fis hmit, Eq. (fin)) reads ^ 
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where in the first line u^fa^^t) is the reduced fioating 
frequency and in the second line we have used that 

Afloat ^ ^Rcs ^ f^S' 

Rescaling the lifetime by typical numbers, we find 
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Notice that, in this case, 6t ^ T^bg, and thus, h cannot 
be treated as a constant. Notice also that 5t scales with 
the —7/3 power of the fioating frequency. This means 
that for even lower floating frequencies, the timescale can 
easily exceed the Hubble time. We find that 5t > Hq^, 
where Hq is the value of the Hubble expansion parame- 
ter today, when 17^^, < 1.7 x IO^'^Mm^hj or equivalently 
/gw ^ 10"'' Hz (IO^Mq/Mmbh), which coincides with 
the low-frequency edge of the LISA sensitivity band. 

Let us now calculate the gravitational wave phase evo- 
lution as an EMRI traverses a floating resonance. We 
apply the second algorithm of Sec. I VII Bl that is bet- 
ter suited to EMRIs, namely that that is derived from 
Eq. dMD: 



ran J 



(116) 



The quantity r can be decomposed into two pieces: (i) 
one that is valid between the initial radius ri„i and the 
radius at which the EMRI begins to float rn^^t; (ii) and 
one that is valid between the floating radius and the final 
radius r^,,. When there is no floating, f — fcR? where the 
latter is given in Eq. (|83p . When there is floating, r 
can be computed considering r = ^(Mmbh, ^Zmbh, ^^orb) in 
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Eq. (l56l) and applying the same chain rule as in Eq. (jlOSp , 

namely 



dr 
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where again a dot denotes a total derivative and in the 
last line we have used Eqs. pTU)) . pTTj) and pT^ . 
The integral in Eq. pi6p can then be solved to obtain 
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The floating radius, the initial radius and the final ra- 
dius can all be written in terms of the time it takes the 
EMRI to traverse a certain distance. When there is no 
floating, we have the standard relation of Eq. (|55|). which 
we can rewrite as 
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where Ti^iji^^t is the time it takes the EMRI to go from 
Tini to Tfioj^t. When there is floating, the relation between 
radius and time changes because f changes: 
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where Tfi„.,t ,:i„ is the time it takes the EMRI to 
go from rji„.,t to r^^. Noting that rfi„^t/MMBH = 
(MMBH^^fioat)"^/^ = (MmbhMs)~^^^i we cau rewrite the 
above expression as 



96 ly"^ 
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The floating time provides information about the val- 
ues of fig for which our treatment here is valid. Since the 
floating time cannot exceed the observation time, and we 
are here setting the latter to one year, then 
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where in the third line we rescaled by typical values. 
Larger values of fis would lead to floating that would ex- 
ceed one-year. Notice that the bound above is perfectly 
consistent with our initial assumption of /^sMmbh <C 1. 
On the other hand, for an EMRI observation to present 
floating resonance modifications, the EMRI had to evolve 
through frequencies inside which there is a fioating reso- 
nance, namely 
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where (/jow'/high) ^^'^ ^^'^ lowest and highest gravita- 
tional wave frequencies of the observed EMRI. We see 
then that there exists a region in fx^ space for which a 
resonance would fall inside of the gravitational wave fre- 
quencies sampled by the EMRI, while at the same time 
leading to floating for less than one year. 

With this at hand, we can rewrite the gravitational 
wave phase of Eq. (|118p in a simpler form. To do so, we 

note that T^^^ ~ T^^^ijio^t + 7ii„.,t_fi„, so that then 
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The dephasing is then nothing but the difference be- tween Eq. (|126p and Eq. ()86p with Tfi„,i„i = T^-bs- Such a 
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FIG. 9: Dephasing in radians for Systems F (top) and R 
(bottom) as a function of the floating time Tfioat.tm, which 
can be mapped to Hs via Eq. p21[) . 



gravitational wave dephasing depends on the usual EMRI 
parameters and also on the floating time Tno^t. Notice 
that as the floating time goes to zero, Eq. (|126p reduces 
exactly the GR gravitational wave phase evolution, and 
thus, the dephasing vanishes. Figure IHl plots this dephas- 
ing as a function of the floating time rfi„.,t for Systems 
F and R. Observe how steep the dephasing is with the 
floating time, reaching values larger than 10^ rads in a 
single hour of floating evolution. 

Such large dephasings would force a template-based, 
matched-filtering search for EMRIs that uses GR tem- 
plates to simply miss the signal all together. The only 
way this dephasing could become small is if fifio^t is very 
close or higher than the highest EMRI frequency ob- 
served: 
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where /cut. high is the maximum frequency resolvable by 
the instrument. If finoat < ■"'/hi™!! then the floating 
resonance would prevent detection of the EMRI in the 
first place. Notice that if Hfl^., < ^fg^, where f^^ ^ 
10~^ Hz is LISA'S low frequency cut-off. then the floating 
would last a time comparable to the Hubble time, and 
thus, the EMRI would never enter the LISA band in the 
first place. 

Given an EMRI observation that is consistent with 
GR, we can then automatically place a constraint on 
the size of /is- This is because if a detection is consis- 
tent with GR, then the floating resonance had to occur 
at a frequency which is extremely close, or higher than. 
Since fin^^t Q.^^, ^ ^s, this implies the con- 
straint /Zs > ■"'/high ^"^^ ^ constrained by the floating 
orbit condition of Eq. ([68)) . 

These results predict strong projected constraints on 
the {^s,ot) phase space shown in Fig. [TJ For example, 
the detection of a System R EMRI would rule out the 
entire region delimited by the curve ABC in Fig. [T] (and 
extending further above and on the left of the plot axis). 



This is an improvement over Solar System constraints of 
over ten orders of magnitude. 



VIII. CONCLUSION AND DISCUSSION 

We have investigated whether scalar-tensor theories 
can be constrained with future EMRI gravitational wave 
observations. The gravitational and scalar emission from 
a EMRI system in a generic (albeit analytic) scalar- 
tensor theory ultimately depends on only two parame- 
ters, (a, /is). Brans-Dicke theory is recovered in a par- 
ticular case when /.t^ = 0. In such theories, scalar dipole 
emission can activate and modify the radiation-reaction 
force via corrections to the energy flux carried away from 
the binary. 

We flrst studied massless scalar tensor theories (/j,s = 
case) and numerically calculated the energy flux carried 
to spatial infinity and into the supermassive black hole 
horizon by the scalar field perturbation. We obtained 
this flux in the Teukolsky formalism, by solving certain 
master equations for a set of given quasi-circular, equa- 
torial geodesies. We then extended the 2.5PN expansion 
of this flux to 3.5PN order through fltting coefficients. 
We also calculated the Fourier transform of the response 
function in the stationary-phase approximation beyond 
Newtonian order. We found that the Brans-Dicke cor- 
rected Fourier phase in the PN expansion can easily be 
mapped to the ppE framework of p^ . [23j and we derived 
the mapping. 

Once the massless scalar energy fiux had been cal- 
culated, we used it in a modified EOB-EMRI frame- 
work to estimate the dephasing induced by the mass- 
less scalar-tensor terms. We found that EMRIs do not 
lead to larger scalar-tensor corrections, but instead, the 
projected bounds on a (or Wbd) are slightly worse than 
those derived from comparable-mass, early inspiral stud- 
ies [1, [3, . In EMRIs, the small companion spends 
all of the observation time very close to the supermas- 
sive black hole. Scalar-tensor theories, however, are not 
necessarily strong-field modifications to GR, but rather 
weak-field ones. That is, the main scalar-tensor modi- 
fication to the orbit is through the introduction of pre- 
Newtonian corrections to the radiation-reaction force due 
to scalar dipole emission, which is dominant for orbits 
with large separations. 

We then studied massive scalar tensor theories (/Xg ^ 
case) and confirmed that certain resonances in the scalar 
horizon flux arise, leading to floating orbits, and thus, 
drastic modifications in the gravitational wave phase evo- 
lution. During floating, the small compact object con- 
tinues to inspiral, but at a rate much slower than the 
standard GR one. The floating inspiral is driven by the 
shedding of mass and spin angular momentum of the su- 
permassive black hole. The timescale for floating to end 
depends on the frequency at which floating occurs, but it 
can easily be greater than 10^ years. In fact, if the EMRI 
separation is greater than ^ 326Mmbh, this timescale ex- 
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ceeds the Hubble time. 

Such large time scales suggest an easy constraint on 
the parameters of massive scalar-tensor theories. That 
is, if these resonances existed, then wc should have found 
many systems orbiting at resonant radii, where material 
that would otherwise be swallowed by the supermassive 
black hole would accumulate. The prototypical system 
that comes to mind are accretion disks, where a super- 
massive black hole accretes material that is driven in- 
wards by viscosity (scalar dipole emission is suppressed 
in neutron star/pulsar binaries due to similar sensitivi- 
ties). Focusing on the orbits of accretion disk material 
in the absence of an EMRI, scalar flux resonances could 
balance the viscous force that the disk would otherwise 
experience, leading to a stagnation of orbits and a build 
up of material at a given radius (see eg. [99l . llOd l). Such 
a build up would translate into a sharp feature in the 
photon energy flux observed, a featu re that in fact is not 
currently needed to fit observations [I03l4l07 |. 

One could argue that this already rules out a certain 
area of parameter space of massive scalar-tensor theories, 
but it is not clear whether such tests can be cleanly ap- 
plied. First, a continuous ring-like, configuration of par- 
ticles in a circular orbit does not lead to superradiance, 
as only the m = modes are excited. More work is nec- 
essary to understand what happens to high (/,m) modes 
that could lead to superradiance when the ring-like con- 
figuration is not continuous, with gaps of the size of the 
orbiting particles. Second, the current observational fits 
of the accretion disk energy spectrum are somewhat sus- 
ceptible to uncertainties in accretion disk phy sics, requir- 
ing the inclusion of non-thermal features |105| and devia- 



tions from the Novikov-Thorne disk model jl07| . There- 
fore, it is not clear whether accretion disks would feel such 
floating resonances, and if so, whether accretion disk ob- 
servations could help us constrain massive scalar-tensor 
theories in a clean manner. 

EMRIs, on the other hand, are much cleaner systems, 
and thus, their observation could be used to place con- 
straints on massive scalar-tensor theories. For a reso- 
nance to affect an EMRI gravitational wave, the former 
must occur at a radius small enough that an orbiting 
EMRI can cross it. In the range of scalar masses in which 
this occurs, the floating orbit introduces an enormous de- 
phasing in the gravitational waveform relative to the GR 
expectation. Such a large dephasing would prevent the 
detection of such modified EMRI waves with LISA or any 
other future detector if one match-filters with GR tem- 
plates. Therefore, given a future EMRI detection that 
is consistent with GR, one could rule out a large area 
of scalar-tensor parameter space. Such projected con- 
straints are several orders of magnitude more stringent 
than current Solar System ones. 

Floating EMRIs emit almost monochromatic gravi- 
tational waves, and thus, one might worry that such 
waves could be detected as a continuous source. These 
waves are almost monochromatic because the radiation- 
reaction force becomes second-order in the mass ratio 



during floating, and since this quantity is much smaller 
than unity, the inspiral evolution is greatly slowed down. 
Continuous algorithms might detect such monochromatic 
waves, but their frequency is much higher than that of 
white-dwarf binaries, the traditional LISA continuous 
sources. In fact, there arc no continuous, monochromatic 
sources that would emit in the high-frequency end of the 
LISA band according to GR. Thus, a continuous source 
detection at a high gravitational wave frequency could 
indicate the presence of floating EMRIs. 

Future work could concentrate on extending the re- 
sults of this paper to relax the condition HsMmsh <C 1. 
Floating orbits can occur near the ISCO of a spinning 
black hole if //sMmbh 0.5. Our results strongly suggest 
that, also in this case, floating orbits would lead to a very 
large dephasing of the gravitational waveform. However, 
a more detailed quantitative analysis, possibly involving 
the numerical fluxes we have here presented, would be 
desirable. One could also investigate whether the non- 
floating resonances have structure and how this modifies 
the waveform observable if the resonant width were large 
enough. 

Another possible avenue for future research would 
be to investigate whether additional constraints can be 
placed on scalar-tensor theories given other astrophys- 
ical observations. For example, if floating orbits exist, 
then the spin of the supermassive black hole would be 
decreased after every EMRI coalescence. Statistically 
speaking, one would then predict a supermassive black 
hole spin distribution that is peaked close to zero. There- 
fore, if a sufficient number of supermassive black holes 
are observed with sufficiently large spins, one would then 
infer that if floating orbits exist, then they had to be hid- 
den inside the horizon of the supermassive black holes. 
This can then be used to place a statistical constr aint on 
the coupling parameters of scalar-tensor theories |l08| . 

Other future work could focus on determining whether 
the floating resonances found here are generic to other 
theories. For example, recently there has been much ef- 
fort to understand certain quadratic gravity theories |14| , 
where the Einstein-Hilbert action is enhanced through 
the addition of terms composed of the product of a scalar 
field and all quadratic curvature invariants. Usually, the 
scalar field in these theories is assumed massless, but if it 
were endowed with a mass, superradiant-induced, float- 
ing resonances could also arise. If so, the same mecha- 
nism described here to constrain the mass of scalar-tensor 
theories could be used to bound the mass of quadratic- 
gravity scalar fields. 
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Appendix A: Flux resonances in the low- frequency, 
low-mass regime 

In this appendix we derive Eq. (|57|) and we generalize 
it for any (Z,m) mode. For simplicity, we set Mmbh = 1 
and all our intermediate derivations are done using di- 
mensionless quantities. At the end of the calculation, our 
results have a simple scaling with Mmbh and it is easy to 
re-insert the necessary factors, as we do in Eq. (|A17[) . 

The procedure is based on solving Eq. by the 

method of matched asymptotic expansions [581 Il09l | . At 
large distance, the wave equation pop is approximately 



0. 



The solution of the equation above with the correct be- 
havior at infinity reads 



where U{a, 6; z) is the hypergeometric function. This so- 
lution can be expanded for small distances as 



AT, 



hnu) 



T[-l - 21] 
V[-l - ijj.'^jkoo] 

"( 2z/loo) 



r[i + 2i] 



-I 



.(Al) 



T[l + 1- iiil/k^] 

Following Detweiler [HI], we define i^^/fcoo = l+l+n+5v. 
Close to the resonances, ik^o ~ —iJ.'^/{l + 1 + n) and 
Siy - 0. Therefore C/(/ + 1 - fil/koo,2l + 2, -2ikoor) ^ 1 
and the wavefunction at the resonance reads 



Imuj 



(A2) 



We have verified this numerically and we find remarkable 
agreement between Eq. (|A2p and the waveform obtained 
by integrating numerically Eq. p9p. all throughout the 
domain. 

On the other hand, close to the horizon, we can ap- 
proximate the wave equation as 



dz 



where Xi, 



.dR 
ziz + 1) — 
dz 



+ [JC^ - l{l + l)z{z + 1)] R , 



,R{z), z={r- r+)/{r+ - r_) 



and K. is defined below Eq. ([57)) . The solution to the 
above equation reads 

R^ {-ly^ {APi{2ilC,l + 2z) + BQi{2iJC,l + 2z)) , 

where Pi and Qi are Legendre polynomials and Legendre 
functions of second kind, respectively. At large distances, 
the near-region wave-function reads 
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Bi-K 3 + e 
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Thus, by matching the two wavefunctions ([Al[) and ([A3 
we can extract the coefficients A,B, which read 
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Now that the solutions have been asymptoticahy 
matched, let us compute the Wronskian. Since the Wron- 
skian of two independent solutions, as defined below 
Eq. (|33|), is constant throughout the space, we can eval- 
uate it at an arbitrary point. In particular, in the near- 
horizon region it reads 



W 



2ifc+Bout , 



(A8) 



where we have used the definition below Eq. P5)) and the 
asymptotic behavior of the homogeneous solutions close 
to the horizon, 



X + ~ p" 



Boute 



ik-^- r« 



(A9) 
(AlO) 



Hence, the Wronskian is simply related to the coefficient 
of the outgoing wave of in the near-horizon re- 

gion. Using Eq. (|A3|) and the coefficients in Eq. (|A6|) 
and (|A7|) ,the solution Ximuj 
becomes 



■^hriLo close to the horizon 



Imuj 



r[i - 2iic] 



A 



2 l-e^^^ 



r.2 I (,2 



iBTi"^ [(1 - tanh"^[27r/C]) sec'^l-K + 2i7r/C]] 

^ 2r[-/ - 2i/c]r[i + 1 - 2iK]v[i + 2iK] ■ 

(All) 

The first term in this expansion is in fact an outgoing 
wave because, close to the horizon. 



(r - r+)" 



so that 



^+'-'{r+ - r^y^ . The refore, Eq. (|jn|) 
has exactly the same form of Eq. (|A10p and one can eas- 
ily extract the coefficient i?out- With this at hand, the 
Wronskian (jA8[) is simply 



r[i - 2i/c] 



X A + 



ittB 1 + e 
~2~ 1 - e*'^" 



and using Eqs. (jA6p and (|A7|) wc find 



-l+iK 



w 
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^ AiV[-~l - 21 - n - 5v] 

which is valid for any {l,m,n). Focusing on the funda- 
mental mode, n = only, at resonance 5v ^ Q and W 
takes the form 



ik+^rl + al^j:[l + l]r[l + I - 2iK] 
" (r+ -r_)-'-^'Cr[2Z + l]r[l-2i/C] ' ^^^^^ 

For / = 1 = m, this simply reduces to 



14^B„= - fc+ ( /C + - ) ^rl + alUr+ - r-f^''^ . (A13) 



Finally, we can now estimate the peak flux coming out 
of the horizon at the resonant frequencies and in the small 
lis limit. Using Eqs. (gH) and (02]), it reads 



peak 



(AM) 



Using ((A2|) we find 



peak 



(A15) 

On the other hand at large distances, from Eq. (|55p . 
'^RcB ^ ^ Ms J s-i^d the expression above simplifies 
to " 

2i-3/2 2M^mo./(i-Hl) 

E;,^ ^ a2^^fc+ |g;,^,(7r/2)p , 

(A16) 

where | Wr^s P is given by Eq. ()A12p . The final expression, 
restoring the necessary factors of A/mbh, reads 



• m/C 7r|5/™(^/2)|2r[l + 2lf (l - ql^^) ' r^/M^^H + S^mbh^WM;^ - 3rH„jA/MBH / r^.., ' 



peak 



4'r[i + ;]2r[i + 1 - 2z/c]r[i + 1 + 2z/c]fc+r+ (^^^^^ + (r,.jM„,H)3/2)' (9Lh + ^L/mSbh) v^-^-- 



We have checked this expression against numerical solu- space for ^sMmbh «C 1. For I = m = 1, Eq. (jA17p re- 
tion and found very good agreement in all of parameter duces to Eq. (|57p in the main text. 
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